This script reads in both PBMC and Liver Seurat objects, and
generates select figures used in the manuscript.
Visualization of metadata
UMAP with no cluster numbers
plot <- DimPlot(sobj, label = FALSE) & NoLegend()
plot
pdf(paste("./figures/", tissue, "/", tissue, "_UMAP_clusters_noLabels.pdf", sep = ""))
plot
dev.off()
png
2

UMAP with cluster numbers
plot <- DimPlot(sobj, label = TRUE)
plot
pdf(paste("./figures/", tissue, "/", tissue, "_UMAP_clusters_labels.pdf", sep = ""))
plot
dev.off()
png
2

Map with SCINA-generated cell-type labels
DimPlot(sobj, group.by = "scina_labels_refined", label = TRUE) & NoLegend()

Map with general cell-type labels for paper
plot <- DimPlot(sobj, group.by = "general_cell_labels",
label = TRUE, repel = TRUE) +
ggtitle(NULL)
plot
pdf(paste("./figures/", tissue, "/", tissue, "_UMAP_general_cell_labels.pdf", sep = ""),
height = 8,
width = 12)
plot
dev.off()
png
2

Map grouping by general cell type labels but with no labels on
plot
plot <- DimPlot(sobj, group.by = "general_cell_labels") +
ggtitle(NULL) &
NoLegend()
plot
pdf(paste("./figures/", tissue, "/", tissue, "_UMAP_general_cell_labels_noLabels.pdf",
sep = ""))
plot
dev.off()
png
2

Map with original identities
plot <- DimPlot(sobj, group.by = "orig.ident",
cols = carto_pal(length(levels(as.factor(sobj$orig.ident))), "Safe")) +
ggtitle(NULL)
plot
pdf(paste("./figures/", tissue, "/", tissue, "_UMAP_orig_idents.pdf",
sep = ""),
height = 5, width = 7)
plot
dev.off()
png
2

Barplot with original identities on a cluster-level grouping:
# Meta data to plot:
df <- sobj@meta.data
# Check what column the cluster identities are in
col <- which(colnames(df) == res)
# Order clusters
df[,col] <- factor(Idents(sobj),
levels = c(sort(as.numeric(levels(Idents(sobj))))))
# Basic plot of clusters by replicate
ggplot(df, aes(x = get(res), fill = orig.ident)) +
geom_bar() +
theme(axis.text = element_text(size = 7))

# Plot as proportion or percentage of cluster
ggplot(df, aes(x = get(res), fill = orig.ident)) +
geom_bar(position = "fill") +
theme(axis.text = element_text(size = 7))

Barplot with original identities on a grouped by general cell
labels:
df <- sobj@meta.data
plot1 <- ggplot(df, aes(x = general_cell_labels, fill = orig.ident)) +
geom_bar() +
scale_fill_carto_d(name = NULL, palette = "Safe") +
theme_bw() +
theme(axis.text = element_text(size = 8),
axis.text.x = element_text(angle = 90, vjust = 0.5, hjust=1),
axis.title.x = element_blank()) +
ylab("Number of cells")
plot1
pdf(paste("./figures/", tissue, "/", tissue, "_barplot_orig_ident_counts.pdf",
sep = ""))
plot1
dev.off()
png
2

# Plot as proportion or percentage of cluster
plot2 <- ggplot(df, aes(x = general_cell_labels, fill = orig.ident)) +
geom_bar(position = "fill") +
scale_fill_carto_d(name = NULL, palette = "Safe") +
theme_bw() +
theme(axis.text = element_text(size = 8),
axis.text.x = element_text(angle = 90, vjust = 0.5, hjust=1),
axis.title.x = element_blank()) +
ylab("Proportion of cells")
plot2
pdf(paste("./figures/", tissue, "/", tissue, "_barplot_orig_ident_proportions.pdf",
sep = ""))
plot2
dev.off()
png
2

Heatmaps
Calculate markers for general cell labels then reset resolution
Idents(sobj) <- "general_cell_labels"
sobj_markers <- RunPrestoAll(sobj,
only.pos = TRUE,
min.pct = 0.25,
logfc.threshold = 0.25)
Calculating cluster B cells
Calculating cluster CD3+/LEF1+/TCF7+ memory T cells
Calculating cluster CD14-/CD16+ Non-classical monocytes
Calculating cluster Erythrocytes
Calculating cluster CD3+/GATA3+/MAF+ TH2-like T cells 1
Calculating cluster CD14-/CD16- Monocytes
Calculating cluster CD3+/CD8A+/NKG7+ NK-like cells
Calculating cluster CD3+/GATA3+/MAF+ TH2-like T cells 2
Calculating cluster CD14+/CD16+ Monocytes
Calculating cluster Antibody-secreting B cells
Calculating cluster Unknown
Calculating cluster Mast cells
Calculating cluster CD3+/MKI67+/TOP2A+ Proliferating T cells
Calculating cluster Hematopoietic stem cells
Calculating cluster Megakaryocytes
sobj_markers %>%
group_by(cluster) %>%
slice_max(n = 2, order_by = avg_log2FC)
Idents(sobj) <- res
Save markers
groups = "general_cell_labels"
write.table(sobj_markers,
file = paste("./figures/", tissue, "/",
tissue, "_markers_", groups, ".tsv",
sep = ""),
quote = FALSE,
sep = "\t",
row.names = FALSE,
col.names = TRUE)
Generate heatmap with top 5 markers grouping by general cell
types
# Remove mikado genes from marker list
sobj_markers <- sobj_markers[grep("mikado", rownames(sobj_markers), invert = TRUE),]
sobj_markers %>%
group_by(cluster) %>%
top_n(n = 5, wt = avg_log2FC) -> top
# If liver, select fewer cells
if (tissue == "liver") {
cells <- sample(colnames(sobj), size = 30000)
} else if (tissue == "PBMC") {
cells <- colnames(sobj)
}
DoHeatmap(sobj, features = top$gene, group.by = "general_cell_labels", size = 3,
angle = 90, cells = cells) +
NoLegend() +
theme(text = element_text(size = 7))

Make PDF of heatmap
groups <- "general_cell_labels"
pdf(paste("./figures/", tissue, "/", tissue, "_heatmap_", groups, ".pdf", sep = ""),
height = 13,
width = 7)
DoHeatmap(sobj, features = top$gene, group.by = "general_cell_labels", size = 2,
angle = 90, cells = cells) +
NoLegend() +
theme(text = element_text(size = 7))
dev.off()
null device
1
Feature plots
Make specific plots with specific genes. The genes we are interested
in include: PTPRC, CALCRL, NKG7, CD3E, MARCO, LYZ, CD19, MS4A1,
STAB2, ALB, CD4, CD8A, CLEC4G, CD5L, C1QB, ACTA2, VWF, IGLL5, CD68.
Can also plot in italics.
geneCode <- "sct_LYZ-1" # Woodchuck-specific nomenclature for this genome
gene <- "LYZ"
mapType <- "Liver"
FeaturePlot(sobj, features = geneCode) +
ggtitle(paste(gene, "-", mapType, "map"))

FeaturePlot(sobj, features = geneCode) +
ggtitle(bquote(~italic(.(gene))))

Another version of the feature plot that outputs genes in italics
if (tissue == "PBMC") {
geneCodes <- c("sct_PTPRC","sct_NKG7","sct_CD14",
"sct_CD3E","sct_MARCO","sct_LYZ-1",
"sct_CD19","sct_MS4A1","sct_STAB2",
"sct_CD4","sct_CD8A", "sct_XCL1;XCL2",
"sct_CD5L","sct_C1QB", "sct_LEF1",
"sct_ACTA2","sct_VWF","sct_IGLL5-1",
"sct_CD68","sct_FCGR3A;FCGR3B","sct_TOP2A")
genes <- c("PTPRC","NKG7","CD14","CD3E","MARCO","LYZ",
"CD19","MS4A1","STAB2","CD4","CD8A", "XCL1;XCL2",
"CD5L","C1QB","LEF1","ACTA2","VWF","IGLL5","CD68",
"FCGR3A;FCGR3B","TOP2A")
} else if (tissue == "liver") {
geneCodes <- c("sct_Ptprc","sct_CALCRL","sct_NKG7",
"sct_CD3E","sct_MARCO","sct_LYZ-1",
"sct_CD19","sct_MS4A1","sct_STAB2",
"sct_ALB-1","sct_CD4","sct_CD8A",
"sct_CLEC4G","sct_CD5L","sct_C1QB",
"sct_ACTA2","sct_VWF","sct_IGLL5-1",
"sct_CD68","sct_XCL1;XCL2","sct_LEF1",
"sct_RSPO3","sct_MECOM")
genes <- c("PTPRC","CALCRL","NKG7","CD3E","MARCO","LYZ",
"CD19","MS4A1","STAB2","ALB","CD4","CD8A",
"CLEC4G","CD5L","C1QB","ACTA2","VWF","IGLL5","CD68",
"XCL1;XCL2","LEF1","RSPO3","MECOM")
}
for(num in 2:length(geneCodes)) {
plot <- FeaturePlot(sobj,
features = geneCodes[num]) +
ggtitle(bquote(~italic(.(genes[num]))))
print(plot)
pdf(paste("./figures/", tissue, "/", tissue, "_", genes[num], "_UMAP.pdf", sep = ""))
print(plot)
dev.off()
}






















Cell-population-specific markers
Plot marker genes for specific cell populations
Endothelial cells:
geneCodes <- c("STAB2","ITGA1","CD55","LYVE1",
"CD34","VWF","IFITM3;IFITM2;IFITM1","RSPO3",
"MECOM","Mecom","CALCRL","LOC114089654",
"RAMP2","BST2","CLEC4G","CTSV;CTSL",
"STAB1","PLAC9","PECAM1","WNT2",
"MYL6","NR2F2")
genes <- c("STAB2","ITGA1","CD55","LYVE1",
"CD34","VWF","IFITM3","RSPO3",
"MECOM","Mecom","CALCRL","LOC114089654",
"RAMP2","BST2","CLEC4G","CTSV;CTSL",
"STAB1","PLAC9","PECAM1","WNT2",
"MYL6","NR2F2")
for(num in 1:length(geneCodes)) {
plot <- FeaturePlot_scCustom(sobj,
features = geneCodes[num],
colors_use = viridis_light_high) +
ggtitle(bquote(~italic(.(genes[num]))))
print(plot)
pdf(paste("./figures/", tissue, "/", tissue, "_", genes[num], "_UMAP.pdf", sep = ""))
print(plot)
dev.off()
}






















Lymphocytes:







Lymphocytes weird genes:

Mesenchyme:
geneCodes <- c("HHIP","COL1A2","COL3A1","IGFBP7",
"IGFBP3","DCN","COL1A1","SPARC",
"RBP1", "CALCRL")
genes <- c("HHIP","COL1A2","COL3A1","IGFBP7",
"IGFBP3","DCN","COL1A1","SPARC",
"RBP1","CALCRL")
for(num in 1:length(geneCodes)) {
plot <- FeaturePlot_scCustom(sobj,
features = geneCodes[num],
colors_use = viridis_light_high) +
ggtitle(bquote(~italic(.(genes[num]))))
print(plot)
pdf(paste("./figures/", tissue, "/", tissue, "_", genes[num], "_UMAP.pdf", sep = ""))
print(plot)
dev.off()
}










Cholangiocytes:
geneCodes <- c("KRT19","CFTR","EPCAM","SLC4A4")
genes <- c("KRT19","CFTR","EPCAM","SLC4A4")
for(num in 1:length(geneCodes)) {
plot <- FeaturePlot_scCustom(sobj,
features = geneCodes[num],
colors_use = viridis_light_high) +
ggtitle(bquote(~italic(.(genes[num]))))
print(plot)
pdf(paste("./figures/", tissue, "/", tissue, "_", genes[num], "_UMAP.pdf", sep = ""))
print(plot)
dev.off()
}




Look at zonation gene signatures:
cv_genes <- c("FETUB","HMGCS1","CYP2E1","GLUD1;GLUD2","CYP1A2",
"RGN","INMT","COMT","FDPS","SPR-1")
cv_genes
[1] "FETUB" "HMGCS1" "CYP2E1" "GLUD1;GLUD2" "CYP1A2" "RGN"
[7] "INMT" "COMT" "FDPS" "SPR-1"
pp_genes <- c("Saa2;Saa1-1","HAMP","APOA1","APOC2","CRYL1",
"AMY1A;AMY1C;AMY1B;AMY2A;AMY2B","MT-ATP6","UROC1",
"APOA2","MT-CO3")
pp_genes
[1] "Saa2;Saa1-1" "HAMP"
[3] "APOA1" "APOC2"
[5] "CRYL1" "AMY1A;AMY1C;AMY1B;AMY2A;AMY2B"
[7] "MT-ATP6" "UROC1"
[9] "APOA2" "MT-CO3"
Cell-type markers all in one plot
Hepatocytes:
geneCodes <- c("CYP2E1","FETUB","HMGCS1","GLUD1;GLUD2",
"CYP1A2","APOC2","Saa1; Saa2", "HAMP",
"PCK1")
plot <- FeaturePlot_scCustom(sobj,
features = geneCodes,
colors_use = viridis_light_high,
num_columns = 4)
plot
pdf(paste("./figures/", tissue, "/", tissue, "_hepatocyteMarkers_UMAP.pdf", sep = ""),
height = ceiling(length(geneCodes)/4)*6, width = 24, width = 24)
plot
dev.off()
Endothelial cells:
# One set of features
geneCodes <- c("STAB2","ITGA1","LYVE1","IFITM3;IFITM2;IFITM1",
"CALCRL","RAMP2","BST2", "CLEC4G",
"CTSV;CTSL", "STAB1", "VWF", "CD34",
"MECOM","PLAC9","RSPO3","WNT2")
# Second set of features
geneCodes <- c("CLEC4G","STAB1","STAB2","LYVE1",
"CD34","MECOM","VWF","RSPO3",
"WNT2","ACKR1")
plot <- FeaturePlot_scCustom(sobj,
features = geneCodes,
colors_use = viridis_light_high,
num_columns = 4)
plot
pdf(paste("./figures/", tissue, "/", tissue, "_endothelialMarkers_UMAP.pdf", sep = ""),
height = ceiling(length(geneCodes)/4)*6, width = 24, width = 24)
plot
dev.off()
T lymphocytes
geneCodes <- c("CD3D","CD3E","CD4","CD8A",
"CTLA4","IL7R-1","LEF1","EOMES",
"TIGIT","KLRB1","KLRD1","NKG7",
"XCL1;XCL2","TOX","GZMA","GZMK")
geneCodes <- c("CD3D","CD3E","GIMAP1","GIMAP4",
"GIMAP6","GIMAP7","CD4","CD28",
"FOXP1","CTLA4","CD38","XCL1;XCL2",
"IL2RA","NKG7","CD8A","CCL5",
"CD69","CD2","GZMK","BRCA2",
"TOP2A", "STMN1")
plot <- FeaturePlot_scCustom(sobj,
features = geneCodes,
colors_use = viridis_light_high,
num_columns = 4)
plot
pdf(paste("./figures/", tissue, "/", tissue, "_lymphocyteMarkers_UMAP.pdf", sep = ""),
height = ceiling(length(geneCodes)/4)*6, width = 24)
plot
dev.off()
T lymphocytes for PBMCs
geneCodes <- c("CD3D","CD3E","XCL1;XCL2","NKG7",
"LEF1","IL7R-1","GATA3","MAF",
"CCR4","S100A4","CXCR3","IL2RA",
"CD160","CD247","BRCA2","STMN1")
# CCR7 and TCF7 aren't showing up???
plot <- FeaturePlot_scCustom(sobj,
features = geneCodes,
colors_use = viridis_light_high,
num_columns = 4,
slot = "data")
plot
pdf(paste("./figures/", tissue, "/", tissue, "_lymphocyteMarkers_UMAP.pdf", sep = ""),
height = ceiling(length(geneCodes)/4)*6, width = 24)
plot
dev.off()
Mesenchyme:
geneCodes <- c("DCN","COL1A2","COL3A1","HHIP",
"CALCRL","SPARC","RBP1","COL1A1",
"IGFBP3", "IGFBP7")
geneCodes <- c("COL1A2", "COL3A1", "IGFBP7", "IGFBP3",
"DCN", "CALCRL", "COL1A1", "SPARC",
"RBP1", "HHIP","RBP1","LRAT",
"PDE3B","HGF","CNN2","ACTA2")
plot <- FeaturePlot_scCustom(sobj,
features = geneCodes,
colors_use = viridis_light_high,
num_columns = 4)
plot
pdf(paste("./figures/", tissue, "/", tissue, "_mesenchymalMarkers_UMAP.pdf", sep = ""),
height = ceiling(length(geneCodes)/4)*6, width = 24)
plot
dev.off()
Cholangiocytes:
geneCodes <- c("KRT19","SLC4A4","CFTR","EPCAM")
geneCodes <- c("ANXA2", "CST3", "BIRC3", "TESC",
"KRT19", "SOX9", "EPCAM","SLC4A4",
"CFTR")
plot <- FeaturePlot_scCustom(sobj,
features = geneCodes,
colors_use = viridis_light_high,
num_columns = 4)
plot
pdf(paste("./figures/", tissue, "/", tissue, "_cholangiocyteMarkers_UMAP.pdf", sep = ""),
height = ceiling(length(geneCodes)/4)*6, width = 24)
plot
dev.off()
Myeloid cells
geneCodes <- c("TYROBP","CD74","CTSS","VSIG4",
"MARCO", "CD5L","HMOX1","VCAM1",
"ITGAM","FCGR3A;FCGR3B","Fcgr3;Fcgr2b","CTSA",
"IL17RA-1","SIRPA;SIRPB1;SIRPG","BTLA", "CADM1",
"ID2","LY6E","FLT3", "CD14",
"NOTCH2","LYZ-1","S100A8")
plot <- FeaturePlot_scCustom(sobj,
features = geneCodes,
colors_use = viridis_light_high,
num_columns = 4)
plot
pdf(paste("./figures/", tissue, "/", tissue, "_myeloidMarkers_UMAP.pdf", sep = ""),
height = ceiling(length(geneCodes)/4)*6, width = 24)
plot
dev.off()
Correlation plots
Generate heatmaps comparing woodchuck clusters with various
datasets
Setup for all correlations
First, set up whatever woodchuck dataset I am working with by reading
in the ortholog table, choosing whether the orthologs to use are human,
mouse, or woodchuck, and calculating the average expression for each
cluster. The output of this section is the scaled cluster
gene-expression matrix
#groups <- res
groups <- "general_cell_labels"
# Read in ortholog table
geneNameTable <- read.table("~/Dropbox/Zoe/scf_version/make_gtf/orthofinder_sc2/homologene/collectedOrthofinderPairings.tsv",
sep = "\t",
header = TRUE)
woodchuckClusterAverages <- AverageExpression(sobj,
assays = "SCT",
slot = "scale.data",
group.by = groups)
# Scale data
woodchuckClusterAverages$SCT <- na.omit(t(scale(t(as.matrix(woodchuckClusterAverages$SCT)))))
# Grab gene names from Seurat object
uniqueHier <- row.names(woodchuckClusterAverages$SCT)
uniqueHier <- as.data.frame(uniqueHier)
# Bind with geneNameTable to get correct order (notice uniqueHier is on left)
newNames <- dplyr::left_join(uniqueHier, geneNameTable, by = "uniqueHier")
# Get orthologs from either mouse or human
species <- "human"
if (species == "human") {
# If human one-to-one ortho has NA, replace with mikado_final_sc2_stringent_noMito_protein column
# This is to avoid and potential mistakes in recognizing things it shouldn't be recognizing
newNames$speciesOneToOne <- ifelse(is.na(newNames$humanOneToOne), newNames$uniqueHier, newNames$humanOneToOne)
} else if (species == "mouse") {
newNames$speciesOneToOne <- ifelse(is.na(newNames$mouseOneToOne), newNames$uniqueHier, newNames$mouseOneToOne)
} else if (species == "woodchuck") {
newNames$speciesOneToOne <- newNames$uniqueHier
}
# Grab dataframe
woodchuckClusterAverages <- woodchuckClusterAverages$SCT
# Replace names with one-to-one orthologue of particular species
row.names(woodchuckClusterAverages) <- newNames$speciesOneToOne
# Make sure formatted correctly
woodchuckClusterAverages <- as.data.frame(woodchuckClusterAverages)
# Order by gene name
woodchuckClusterAverages <- woodchuckClusterAverages[order(row.names(woodchuckClusterAverages)),]
# Sanity check
head(woodchuckClusterAverages)
Dataset-specific setup
Correlation of woodchuck PBMCs with human 68k PBMC dataset from 10X
Genomics
# Now need to read in 68K PBMC data
humanPBMC <- read.csv("~/Dropbox/Zoe/scf_version/analysis/correlationTests/68K_pbmc_data/68K_enrichedGenes.csv",
header = FALSE)
# Get rid of top row
humanPBMC <- humanPBMC[-1,]
# Separate all cell and myloid cell data
allCells <- select(humanPBMC, V1, V2, V3)
myeloid <- select(humanPBMC, V5, V6, V7)
# Rename and get rid of first row
colnames(allCells) <- c("Cluster", "Gene", "Enrichment")
colnames(myeloid) <- c("Cluster", "Gene", "Enrichment")
# Get rid of top row
allCells <- allCells[-1,]
myeloid <- myeloid[-1,]
# I can then filter the rows and bind the data frames back together by gene name
for (j in 1:10) {
dat <- dplyr::filter(allCells, Cluster == j)
#dat <- get(paste("allCells", j, sep = ""))
# Multiply enrichment values by -1 because the signs are backwards???
dat$Enrichment <- as.numeric(dat$Enrichment) * -1
dat <- dplyr::select(dat, Gene, Enrichment)
colnames(dat) <- c("Gene", paste("Enrichment", j, sep = ""))
if (j == 1) {
allCellsMatrix <- dat
}
else {
allCellsMatrix <- dplyr::full_join(allCellsMatrix, dat, by = "Gene")
}
}
# Make row names gene names
rownames(allCellsMatrix) <- allCellsMatrix$Gene
allCellsMatrix <- dplyr::select(allCellsMatrix, -Gene)
# Make column names cell types
colnames(allCellsMatrix) <- c("Activated CD8+", "Naive CD8+", "Memory and Reg T",
"Naive CD4+", "NK", "CD8+", "B", "Megakaryocytes",
"Monocytes and Dendritic", "B, Dendritic, T")
# Scale across columns
allCellsMatrix <- t(scale(t(allCellsMatrix)))
# Order genes alphabetically by gene name
allCellsMatrix <- allCellsMatrix[order(row.names(allCellsMatrix)),]
# Sanity check
head(allCellsMatrix)
speciesData <- "68KPBMC"
Correlation of woodchuck liver with human liver dataset from
MacParland et al. (2018)
# Find cluster averages of human liver data
load("~/Dropbox/Zoe/scf_version/analysis/correlationTests/HumanLiver.RData")
# Run SCTransform
HumanLiverSeurat <- UpdateSeuratObject(HumanLiverSeurat)
HumanLiverSeurat <- SCTransform(HumanLiverSeurat)
humanClusterAverages <- AverageExpression(HumanLiverSeurat,
assays = "SCT",
slot = "scale.data")
# Replace cluster numbers with names
colnames(humanClusterAverages$SCT) <- c("Hep 1", "Alpha-beta T cells", "Hep 2",
"Inflammatory macs", "Hep 3", "Hep 4",
"Plasma cells", "NK-like cells", "Gamma-delta T cells",
"Non-inflammatory macs", "Periportal LSECs", "Central venous LSECs",
"Portal endothelial cells", "Hep 5", "Hep 6",
"Mature B cells", "Cholangiocytes", "Gamma-delta T cells 2",
"Erythroid cells", "Hepatic stellate cells")
# If only looking at specific clusters
#humanClusterAverages$SCT <- humanClusterAverages$SCT[,c("3","1","15","6","14","5")]
# Otherwise go straight to here:
humanClusterAverages$SCT <- na.omit(t(scale(t(as.matrix(humanClusterAverages$SCT)))))
# Grab gene names
humanGenes <- row.names(humanClusterAverages$SCT)
# Now turn into large dataframe
allCellsMatrix <- as.data.frame(humanClusterAverages$SCT)
# Order by row name
allCellsMatrix <- allCellsMatrix[order(row.names(allCellsMatrix)),]
# Sanity check
head(allCellsMatrix)
speciesData <- "macparland"
Correlation of woodchuck liver with human liver dataset from Aizarani
et al.
# Read in Aizarani dataset
aizarani <- readRDS("~/Dropbox/Zoe/scf_version/analysis/correlationTests/GSE124395_Normalhumanliverdata.RData")
# Read in clusters and label cells
aizaraniClusters <- read.table("~/Dropbox/Zoe/scf_version/analysis/correlationTests/GSE124395_clusterpartition.txt")
# Only keep cells in the cluster object
aizarani <- aizarani[,intersect(colnames(aizarani),row.names(aizaraniClusters))]
# Create Seurat object
aizarani <- CreateSeuratObject(counts = aizarani)
# Run SCTransform
aizarani <- SCTransform(aizarani)
# Add cluster IDs
Idents(aizarani) <- aizaraniClusters$sct.cpart
# Get cluster averages
aizaraniAverages <- AverageExpression(aizarani,
assays = "SCT",
slot = "scale.data")
aizaraniAverages$SCT <- na.omit(t(scale(t(as.matrix(aizaraniAverages$SCT)))))
# Grab gene names
aizaraniGenes <- row.names(aizaraniAverages$SCT)
# Now turn into large dataframe
allCellsMatrix <- as.data.frame(aizaraniAverages$SCT)
# Order by column
allCellsMatrix <- allCellsMatrix[,as.character(sort(as.numeric(colnames(allCellsMatrix))))]
# Rename columns to be more meaningful (not totally confident I got all correct)
colnames(allCellsMatrix) <- c('NK, NKT and T cells (1)',
'Kupffer cells (2)',
'NK, NKT and T cells (3)',
'EPCAM+ cells and cholangiocytes (4)',
'NK, NKT and T cells (5)',
'Kupffer cells (6)',
'EPCAM+ cells and cholangiocytes (7)',
'B cells (8)',
'Liver sinusoidal endothelial cells (9)',
'Macrovascular endothelial cells (10)',
'Hepatocytes (11)',
'NK, NKT and T cells (12)',
'Liver sinusoidal endothelial cells (13)',
'Hepatocytes (14)',
'Other endothelial cells (15)',
'Other (16)',
'Hepatocytes (17)',
'NK, NKT and T cells (18)',
'NK, NKT and T cells (19)',
'Liver sinusoidal endothelial cells (20)',
'Macrovascular endothelial cells (21)',
'B cells (22)',
'Kupffer cells (23)',
'EPCAM+ cells and cholangiocytes (24)',
'Kupffer cells (25)',
'Other endothelial cells (26)',
'Other (27)',
'NK, NKT and T cells (28)',
'Macrovascular endothelial cells (29)',
'Hepatocytes (30)',
'Kupffer cells (31)',
'Macrovascular endothelial cells (32)',
'Stellate cells and myofibroblasts (33)',
'B cells (34)',
'Other endothelial cells (35)',
'Other (36)',
'Other (37)',
'B cells (38)',
'EPCAM+ cells and cholangiocytes (39)')
# Order by row name
allCellsMatrix <- allCellsMatrix[order(row.names(allCellsMatrix)),]
# Sanity check
head(allCellsMatrix)
speciesData <- "aizarani"
Correlation of woodchuck liver with human liver dataset from Andrews
et al. (2022)
# Load dataset
humanNuc <- readRDS("~/Dropbox/Zoe/scf_version/analysis/correlationTests/single_nuc_20_human_map.rds")
# Isolate single-nuc only
humanNuc <- UpdateSeuratObject(humanNuc)
humanNuc <- subset(humanNuc, subset = assay_type == "single_nuc")
# Run SCTransform
humanNuc <- SCTransform(humanNuc)
# Change idents
Idents(humanNuc) <- humanNuc@meta.data$sub_annotation
# Get cluster averages
humanNucClustAverages <- AverageExpression(humanNuc,
assays = "SCT",
slot = "scale.data")
humanNucClustAverages$SCT <- na.omit(t(scale(t(as.matrix(humanNucClustAverages$SCT)))))
# Grab gene names
nucGenes <- row.names(humanNucClustAverages$SCT)
# Now turn into large dataframe
allCellsMatrix <- as.data.frame(humanNucClustAverages$SCT)
# Order by row name
allCellsMatrix <- allCellsMatrix[order(row.names(allCellsMatrix)),]
# Sanity check
head(allCellsMatrix)
speciesData <- "andrews"
Correlation of woodchuck liver with woodchuck PBMCs. For this
correlation, read in the woodchuck liver dataset at the beginning of
this script and then read in the woodchuck PBMCs below
groups <- "general_cell_labels"
# Start with liver and read in woodchuck PBMCs again
load("~/Dropbox/Zoe/scf_version/analysis/healthy_sc/seurat_objects/no_dropletQC/integrated_PBMC_cca_kanchor5_scClustViz.RData")
Idents(scSeurat) <- "integrated_snn_res.0.6"
# Find cluster averages
pbmcClusterAverages <- AverageExpression(scSeurat,
assays = "SCT",
slot = "scale.data",
group.by = groups)
pbmcClusterAverages <- as.data.frame(na.omit(t(scale(t(as.matrix(pbmcClusterAverages$SCT))))))
# Order by row name
allCellsMatrix <- pbmcClusterAverages[order(row.names(pbmcClusterAverages)),]
speciesData <- "PBMC"
Output plots for all correlations
# Now find intersecting genes
matches <- intersect(row.names(allCellsMatrix),
row.names(woodchuckClusterAverages))
# Look at how many genes matched
length(matches)
# Make new matrices with only matching gene names
toCor <- allCellsMatrix[matches,]
woodchuckAveragesCor <- woodchuckClusterAverages[matches,]
# Do Pearson
pearVal <- cor(toCor, woodchuckAveragesCor, method = "pearson")
heatmap(pearVal)
#main = paste("Pearson correlation of", speciesData, "vs", woodchuckData),
#xlab = woodchuckData,
#ylab = speciesData)
#margins = c(6,11))
#Rowv = NA,
#Colv = NA)
pdf(paste("./figures/", tissue, "/", tissue, "_", speciesData, "_PearsonCor.pdf", sep = ""),
height = 15, width = 13)
heatmap(pearVal, margins = c(13,13))
dev.off()
# Do Spearman
spearVal <- cor(toCor, woodchuckAveragesCor, method = "spearman")
heatmap(spearVal)
#main = paste("Spearman correlation of", speciesData, "vs", woodchuckData),
#xlab = woodchuckData,
#ylab = speciesData)
#margins = c(6,11))
#Rowv = NA,
#Colv = NA)
pdf(paste("./figures/", tissue, "/", tissue, "_", speciesData, "_SpearmanCor.pdf", sep = ""),
height = 15, width = 13)
heatmap(spearVal, margins = c(13,13))
dev.off()
LS0tCnRpdGxlOiAiUiBOb3RlYm9vayIKb3V0cHV0OiBodG1sX25vdGVib29rCi0tLQoKVGhpcyBzY3JpcHQgcmVhZHMgaW4gYm90aCBQQk1DIGFuZCBMaXZlciBTZXVyYXQgb2JqZWN0cywgYW5kIGdlbmVyYXRlcyBzZWxlY3QgZmlndXJlcyB1c2VkIGluIHRoZSBtYW51c2NyaXB0LgoKIyMgUHJlcGFyZSBvYmplY3RzCgpMb2FkIGxpYnJhcmllcwoKYGBge3J9CmxpYnJhcnkoU2V1cmF0KQpsaWJyYXJ5KHNjQ2x1c3RWaXopCmxpYnJhcnkoZ2dwbG90MikKbGlicmFyeShkcGx5cikKbGlicmFyeShyY2FydG9jb2xvcikKbGlicmFyeShTZXVyYXRXcmFwcGVycykKbGlicmFyeShzY0N1c3RvbWl6ZSkKYGBgCgpSZWFkIGluIGxpdmVyIG1hcAoKYGBge3J9CnNvYmogPC0gcmVhZFJEUygifi9Ecm9wYm94L1pvZS9zY2ZfdmVyc2lvbi9hbmFseXNpcy9oZWFsdGh5X3NjL3NldXJhdF9vYmplY3RzL2Ryb3BsZXRRQ19maWx0ZXJlZC9hbGxJbnRlZ3JhdGVkX2NjYV9rYW5jaG9yNV9ub0Jpb3BzeUhlcHNfZHJvcGxldFFDRmlsdGVyZWQuUkRTIikKcmVzIDwtICJpbnRlZ3JhdGVkX3Nubl9yZXMuMS40IgpJZGVudHMoc29iaikgPC0gcmVzCnRpc3N1ZSA8LSAibGl2ZXIiCmBgYAoKUmVhZCBpbiBQQk1DIG1hcAoKYGBge3J9CmxvYWQoIn4vRHJvcGJveC9ab2Uvc2NmX3ZlcnNpb24vYW5hbHlzaXMvaGVhbHRoeV9zYy9zZXVyYXRfb2JqZWN0cy9ub19kcm9wbGV0UUMvaW50ZWdyYXRlZF9QQk1DX2NjYV9rYW5jaG9yNV9zY0NsdXN0Vml6LlJEYXRhIikKc29iaiA8LSBzY1NldXJhdApyZXMgPC0gImludGVncmF0ZWRfc25uX3Jlcy4wLjYiCklkZW50cyhzb2JqKSA8LSByZXMKdGlzc3VlIDwtICJQQk1DIgpgYGAKCiMjIFZpc3VhbGl6YXRpb24gb2YgbWV0YWRhdGEKClVNQVAgd2l0aCBubyBjbHVzdGVyIG51bWJlcnMKCmBgYHtyfQpwbG90IDwtIERpbVBsb3Qoc29iaiwgbGFiZWwgPSBGQUxTRSkgJiBOb0xlZ2VuZCgpCnBsb3QKcGRmKHBhc3RlKCIuL2ZpZ3VyZXMvIiwgdGlzc3VlLCAiLyIsIHRpc3N1ZSwgIl9VTUFQX2NsdXN0ZXJzX25vTGFiZWxzLnBkZiIsIHNlcCA9ICIiKSkKcGxvdApkZXYub2ZmKCkKYGBgCgpVTUFQIHdpdGggY2x1c3RlciBudW1iZXJzCgpgYGB7cn0KcGxvdCA8LSBEaW1QbG90KHNvYmosIGxhYmVsID0gVFJVRSkKcGxvdApwZGYocGFzdGUoIi4vZmlndXJlcy8iLCB0aXNzdWUsICIvIiwgdGlzc3VlLCAiX1VNQVBfY2x1c3RlcnNfbGFiZWxzLnBkZiIsIHNlcCA9ICIiKSkKcGxvdApkZXYub2ZmKCkKYGBgCgpNYXAgd2l0aCBTQ0lOQS1nZW5lcmF0ZWQgY2VsbC10eXBlIGxhYmVscwoKYGBge3J9CkRpbVBsb3Qoc29iaiwgZ3JvdXAuYnkgPSAic2NpbmFfbGFiZWxzX3JlZmluZWQiLCBsYWJlbCA9IFRSVUUpICYgTm9MZWdlbmQoKQpgYGAKCk1hcCB3aXRoIGdlbmVyYWwgY2VsbC10eXBlIGxhYmVscyBmb3IgcGFwZXIKCmBgYHtyfQpwbG90IDwtIERpbVBsb3Qoc29iaiwgZ3JvdXAuYnkgPSAiZ2VuZXJhbF9jZWxsX2xhYmVscyIsCiAgICAgICAgbGFiZWwgPSBUUlVFLCByZXBlbCA9IFRSVUUpICsKICBnZ3RpdGxlKE5VTEwpCnBsb3QKcGRmKHBhc3RlKCIuL2ZpZ3VyZXMvIiwgdGlzc3VlLCAiLyIsIHRpc3N1ZSwgIl9VTUFQX2dlbmVyYWxfY2VsbF9sYWJlbHMucGRmIiwgc2VwID0gIiIpLAogICAgaGVpZ2h0ID0gOCwKICAgIHdpZHRoID0gMTIpCnBsb3QKZGV2Lm9mZigpCmBgYAoKTWFwIGdyb3VwaW5nIGJ5IGdlbmVyYWwgY2VsbCB0eXBlIGxhYmVscyBidXQgd2l0aCBubyBsYWJlbHMgb24gcGxvdAoKYGBge3J9CnBsb3QgPC0gRGltUGxvdChzb2JqLCBncm91cC5ieSA9ICJnZW5lcmFsX2NlbGxfbGFiZWxzIikgKwogIGdndGl0bGUoTlVMTCkgJgogIE5vTGVnZW5kKCkKcGxvdApwZGYocGFzdGUoIi4vZmlndXJlcy8iLCB0aXNzdWUsICIvIiwgdGlzc3VlLCAiX1VNQVBfZ2VuZXJhbF9jZWxsX2xhYmVsc19ub0xhYmVscy5wZGYiLCAKICAgICAgICAgIHNlcCA9ICIiKSkKcGxvdApkZXYub2ZmKCkKYGBgCgpNYXAgd2l0aCBvcmlnaW5hbCBpZGVudGl0aWVzCgpgYGB7cn0KcGxvdCA8LSBEaW1QbG90KHNvYmosIGdyb3VwLmJ5ID0gIm9yaWcuaWRlbnQiLAogICAgICAgIGNvbHMgPSBjYXJ0b19wYWwobGVuZ3RoKGxldmVscyhhcy5mYWN0b3Ioc29iaiRvcmlnLmlkZW50KSkpLCAiU2FmZSIpKSArCiAgZ2d0aXRsZShOVUxMKQpwbG90CnBkZihwYXN0ZSgiLi9maWd1cmVzLyIsIHRpc3N1ZSwgIi8iLCB0aXNzdWUsICJfVU1BUF9vcmlnX2lkZW50cy5wZGYiLCAKICAgICAgICAgIHNlcCA9ICIiKSwKICAgIGhlaWdodCA9IDUsIHdpZHRoID0gNykKcGxvdApkZXYub2ZmKCkKYGBgCgpCYXJwbG90IHdpdGggb3JpZ2luYWwgaWRlbnRpdGllcyBvbiBhIGNsdXN0ZXItbGV2ZWwgZ3JvdXBpbmc6CgpgYGB7cn0KIyBNZXRhIGRhdGEgdG8gcGxvdDoKZGYgPC0gc29iakBtZXRhLmRhdGEKIyBDaGVjayB3aGF0IGNvbHVtbiB0aGUgY2x1c3RlciBpZGVudGl0aWVzIGFyZSBpbgpjb2wgPC0gd2hpY2goY29sbmFtZXMoZGYpID09IHJlcykKIyBPcmRlciBjbHVzdGVycwpkZlssY29sXSA8LSBmYWN0b3IoSWRlbnRzKHNvYmopLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgIGxldmVscyA9IGMoc29ydChhcy5udW1lcmljKGxldmVscyhJZGVudHMoc29iaikpKSkpKQojIEJhc2ljIHBsb3Qgb2YgY2x1c3RlcnMgYnkgcmVwbGljYXRlCmdncGxvdChkZiwgYWVzKHggPSBnZXQocmVzKSwgZmlsbCA9IG9yaWcuaWRlbnQpKSArCiAgZ2VvbV9iYXIoKSArCiAgdGhlbWUoYXhpcy50ZXh0ID0gZWxlbWVudF90ZXh0KHNpemUgPSA3KSkKIyBQbG90IGFzIHByb3BvcnRpb24gb3IgcGVyY2VudGFnZSBvZiBjbHVzdGVyCmdncGxvdChkZiwgYWVzKHggPSBnZXQocmVzKSwgZmlsbCA9IG9yaWcuaWRlbnQpKSArCiAgZ2VvbV9iYXIocG9zaXRpb24gPSAiZmlsbCIpICsKICB0aGVtZShheGlzLnRleHQgPSBlbGVtZW50X3RleHQoc2l6ZSA9IDcpKQpgYGAKCkJhcnBsb3Qgd2l0aCBvcmlnaW5hbCBpZGVudGl0aWVzIG9uIGEgZ3JvdXBlZCBieSBnZW5lcmFsIGNlbGwgbGFiZWxzOgoKYGBge3J9CmRmIDwtIHNvYmpAbWV0YS5kYXRhCnBsb3QxIDwtIGdncGxvdChkZiwgYWVzKHggPSBnZW5lcmFsX2NlbGxfbGFiZWxzLCBmaWxsID0gb3JpZy5pZGVudCkpICsKICBnZW9tX2JhcigpICsKICBzY2FsZV9maWxsX2NhcnRvX2QobmFtZSA9IE5VTEwsIHBhbGV0dGUgPSAiU2FmZSIpICsKICB0aGVtZV9idygpICsKICB0aGVtZShheGlzLnRleHQgPSBlbGVtZW50X3RleHQoc2l6ZSA9IDgpLAogICAgICAgIGF4aXMudGV4dC54ID0gZWxlbWVudF90ZXh0KGFuZ2xlID0gOTAsIHZqdXN0ID0gMC41LCBoanVzdD0xKSwKICAgICAgICBheGlzLnRpdGxlLnggPSBlbGVtZW50X2JsYW5rKCkpICsKICB5bGFiKCJOdW1iZXIgb2YgY2VsbHMiKQpwbG90MQpwZGYocGFzdGUoIi4vZmlndXJlcy8iLCB0aXNzdWUsICIvIiwgdGlzc3VlLCAiX2JhcnBsb3Rfb3JpZ19pZGVudF9jb3VudHMucGRmIiwgCiAgICAgICAgICBzZXAgPSAiIikpCnBsb3QxCmRldi5vZmYoKQojIFBsb3QgYXMgcHJvcG9ydGlvbiBvciBwZXJjZW50YWdlIG9mIGNsdXN0ZXIKcGxvdDIgPC0gZ2dwbG90KGRmLCBhZXMoeCA9IGdlbmVyYWxfY2VsbF9sYWJlbHMsIGZpbGwgPSBvcmlnLmlkZW50KSkgKwogIGdlb21fYmFyKHBvc2l0aW9uID0gImZpbGwiKSArCiAgc2NhbGVfZmlsbF9jYXJ0b19kKG5hbWUgPSBOVUxMLCBwYWxldHRlID0gIlNhZmUiKSArCiAgdGhlbWVfYncoKSArCiAgdGhlbWUoYXhpcy50ZXh0ID0gZWxlbWVudF90ZXh0KHNpemUgPSA4KSwKICAgICAgICBheGlzLnRleHQueCA9IGVsZW1lbnRfdGV4dChhbmdsZSA9IDkwLCB2anVzdCA9IDAuNSwgaGp1c3Q9MSksCiAgICAgICAgYXhpcy50aXRsZS54ID0gZWxlbWVudF9ibGFuaygpKSArCiAgeWxhYigiUHJvcG9ydGlvbiBvZiBjZWxscyIpCnBsb3QyCnBkZihwYXN0ZSgiLi9maWd1cmVzLyIsIHRpc3N1ZSwgIi8iLCB0aXNzdWUsICJfYmFycGxvdF9vcmlnX2lkZW50X3Byb3BvcnRpb25zLnBkZiIsIAogICAgICAgICAgc2VwID0gIiIpKQpwbG90MgpkZXYub2ZmKCkKYGBgCgojIyBEb3RwbG90cwoKR2VuZXJhdGUgZG90cGxvdCB3aXRoIHNwZWNpZmljIG1hcmtlcnMKCmBgYHtyfQpEb3RQbG90KHNvYmosCiAgICAgICAgYXNzYXkgPSAiU0NUIiwKICAgICAgICBmZWF0dXJlcyA9IGMoIlBUUFJDIiwgIkNBTENSTCIsICJOS0c3IiwgIkNEM0UiLCAiTUFSQ08iLCAiTFlaLTEiLCAiQ0QxOSIsICJNUzRBMSIsICJTVEFCMiIpCiAgICAgICAgKSArCiAgZ2d0aXRsZSgiU2VsZWN0IGZlYXR1cmVzIGZvciBsaXZlciBtYXAiKQpgYGAKCiMjIEhlYXRtYXBzCgpDYWxjdWxhdGUgbWFya2VycyBmb3IgZ2VuZXJhbCBjZWxsIGxhYmVscyB0aGVuIHJlc2V0IHJlc29sdXRpb24KCmBgYHtyfQpJZGVudHMoc29iaikgPC0gImdlbmVyYWxfY2VsbF9sYWJlbHMiCnNvYmpfbWFya2VycyA8LSBSdW5QcmVzdG9BbGwoc29iaiwgCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgb25seS5wb3MgPSBUUlVFLCAKICAgICAgICAgICAgICAgICAgICAgICAgICAgICBtaW4ucGN0ID0gMC4yNSwgCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgbG9nZmMudGhyZXNob2xkID0gMC4yNSkKc29ial9tYXJrZXJzICU+JQogICAgZ3JvdXBfYnkoY2x1c3RlcikgJT4lCiAgICBzbGljZV9tYXgobiA9IDIsIG9yZGVyX2J5ID0gYXZnX2xvZzJGQykKSWRlbnRzKHNvYmopIDwtIHJlcwpgYGAKClNhdmUgbWFya2VycwoKYGBge3J9Cmdyb3VwcyA9ICJnZW5lcmFsX2NlbGxfbGFiZWxzIgp3cml0ZS50YWJsZShzb2JqX21hcmtlcnMsCiAgICAgICAgICAgIGZpbGUgPSBwYXN0ZSgiLi9maWd1cmVzLyIsIHRpc3N1ZSwgIi8iLAogICAgICAgICAgICAgICAgICAgICAgICAgdGlzc3VlLCAiX21hcmtlcnNfIiwgZ3JvdXBzLCAiLnRzdiIsIAogICAgICAgICAgICAgICAgICAgICAgICAgc2VwID0gIiIpLAogICAgICAgICAgICBxdW90ZSA9IEZBTFNFLAogICAgICAgICAgICBzZXAgPSAiXHQiLAogICAgICAgICAgICByb3cubmFtZXMgPSBGQUxTRSwKICAgICAgICAgICAgY29sLm5hbWVzID0gVFJVRSkKYGBgCgpHZW5lcmF0ZSBoZWF0bWFwIHdpdGggdG9wIDUgbWFya2VycyBncm91cGluZyBieSBnZW5lcmFsIGNlbGwgdHlwZXMKCmBgYHtyfQojIFJlbW92ZSBtaWthZG8gZ2VuZXMgZnJvbSBtYXJrZXIgbGlzdApzb2JqX21hcmtlcnMgPC0gc29ial9tYXJrZXJzW2dyZXAoIm1pa2FkbyIsIHJvd25hbWVzKHNvYmpfbWFya2VycyksIGludmVydCA9IFRSVUUpLF0Kc29ial9tYXJrZXJzICU+JQogICAgZ3JvdXBfYnkoY2x1c3RlcikgJT4lCiAgICB0b3BfbihuID0gNSwgd3QgPSBhdmdfbG9nMkZDKSAtPiB0b3AKIyBJZiBsaXZlciwgc2VsZWN0IGZld2VyIGNlbGxzCmlmICh0aXNzdWUgPT0gImxpdmVyIikgewogIGNlbGxzIDwtIHNhbXBsZShjb2xuYW1lcyhzb2JqKSwgc2l6ZSA9IDMwMDAwKQp9IGVsc2UgaWYgKHRpc3N1ZSA9PSAiUEJNQyIpIHsKICBjZWxscyA8LSBjb2xuYW1lcyhzb2JqKQp9CkRvSGVhdG1hcChzb2JqLCBmZWF0dXJlcyA9IHRvcCRnZW5lLCBncm91cC5ieSA9ICJnZW5lcmFsX2NlbGxfbGFiZWxzIiwgc2l6ZSA9IDMsIAogICAgICAgICAgYW5nbGUgPSA5MCwgY2VsbHMgPSBjZWxscykgKwogIE5vTGVnZW5kKCkgKwogIHRoZW1lKHRleHQgPSBlbGVtZW50X3RleHQoc2l6ZSA9IDcpKQpgYGAKCk1ha2UgUERGIG9mIGhlYXRtYXAKCmBgYHtyfQpncm91cHMgPC0gImdlbmVyYWxfY2VsbF9sYWJlbHMiCnBkZihwYXN0ZSgiLi9maWd1cmVzLyIsIHRpc3N1ZSwgIi8iLCB0aXNzdWUsICJfaGVhdG1hcF8iLCBncm91cHMsICIucGRmIiwgc2VwID0gIiIpLAogICAgaGVpZ2h0ID0gMTMsCiAgICB3aWR0aCA9IDcpCkRvSGVhdG1hcChzb2JqLCBmZWF0dXJlcyA9IHRvcCRnZW5lLCBncm91cC5ieSA9ICJnZW5lcmFsX2NlbGxfbGFiZWxzIiwgc2l6ZSA9IDIsIAogICAgICAgICAgYW5nbGUgPSA5MCwgY2VsbHMgPSBjZWxscykgKwogIE5vTGVnZW5kKCkgKwogIHRoZW1lKHRleHQgPSBlbGVtZW50X3RleHQoc2l6ZSA9IDcpKQpkZXYub2ZmKCkKYGBgCgoKIyMgRmVhdHVyZSBwbG90cwoKTWFrZSBzcGVjaWZpYyBwbG90cyB3aXRoIHNwZWNpZmljIGdlbmVzLiBUaGUgZ2VuZXMgd2UgYXJlIGludGVyZXN0ZWQgaW4gaW5jbHVkZTogKlBUUFJDLCBDQUxDUkwsIE5LRzcsIENEM0UsIE1BUkNPLCBMWVosIENEMTksIE1TNEExLCBTVEFCMiwgQUxCLCBDRDQsIENEOEEsIENMRUM0RywgQ0Q1TCwgQzFRQiwgQUNUQTIsIFZXRiwgSUdMTDUsIENENjgqLiBDYW4gYWxzbyBwbG90IGluIGl0YWxpY3MuCgpgYGB7cn0KZ2VuZUNvZGUgPC0gInNjdF9MWVotMSIgIyBXb29kY2h1Y2stc3BlY2lmaWMgbm9tZW5jbGF0dXJlIGZvciB0aGlzIGdlbm9tZQpnZW5lIDwtICJMWVoiCm1hcFR5cGUgPC0gIkxpdmVyIgpGZWF0dXJlUGxvdChzb2JqLCBmZWF0dXJlcyA9IGdlbmVDb2RlKSArCiAgZ2d0aXRsZShwYXN0ZShnZW5lLCAiLSIsIG1hcFR5cGUsICJtYXAiKSkKRmVhdHVyZVBsb3Qoc29iaiwgZmVhdHVyZXMgPSBnZW5lQ29kZSkgKwogIGdndGl0bGUoYnF1b3RlKH5pdGFsaWMoLihnZW5lKSkpKQpgYGAKCkFub3RoZXIgdmVyc2lvbiBvZiB0aGUgZmVhdHVyZSBwbG90IHRoYXQgb3V0cHV0cyBnZW5lcyBpbiBpdGFsaWNzCgpgYGB7cn0KaWYgKHRpc3N1ZSA9PSAiUEJNQyIpIHsKIGdlbmVDb2RlcyA8LSBjKCJzY3RfUFRQUkMiLCJzY3RfTktHNyIsInNjdF9DRDE0IiwKICAgICAgICAgICAgICAgInNjdF9DRDNFIiwic2N0X01BUkNPIiwic2N0X0xZWi0xIiwKICAgICAgICAgICAgICAgInNjdF9DRDE5Iiwic2N0X01TNEExIiwic2N0X1NUQUIyIiwKICAgICAgICAgICAgICAgInNjdF9DRDQiLCJzY3RfQ0Q4QSIsICJzY3RfWENMMTtYQ0wyIiwKICAgICAgICAgICAgICAgInNjdF9DRDVMIiwic2N0X0MxUUIiLCAic2N0X0xFRjEiLAogICAgICAgICAgICAgICAic2N0X0FDVEEyIiwic2N0X1ZXRiIsInNjdF9JR0xMNS0xIiwKICAgICAgICAgICAgICAgInNjdF9DRDY4Iiwic2N0X0ZDR1IzQTtGQ0dSM0IiLCJzY3RfVE9QMkEiKQogZ2VuZXMgPC0gYygiUFRQUkMiLCJOS0c3IiwiQ0QxNCIsIkNEM0UiLCJNQVJDTyIsIkxZWiIsCiAgICAgICAgICAgIkNEMTkiLCJNUzRBMSIsIlNUQUIyIiwiQ0Q0IiwiQ0Q4QSIsICJYQ0wxO1hDTDIiLAogICAgICAgICAgICJDRDVMIiwiQzFRQiIsIkxFRjEiLCJBQ1RBMiIsIlZXRiIsIklHTEw1IiwiQ0Q2OCIsCiAgICAgICAgICAgIkZDR1IzQTtGQ0dSM0IiLCJUT1AyQSIpIAp9IGVsc2UgaWYgKHRpc3N1ZSA9PSAibGl2ZXIiKSB7CiAgZ2VuZUNvZGVzIDwtIGMoInNjdF9QdHByYyIsInNjdF9DQUxDUkwiLCJzY3RfTktHNyIsCiAgICAgICAgICAgICAgICAgInNjdF9DRDNFIiwic2N0X01BUkNPIiwic2N0X0xZWi0xIiwKICAgICAgICAgICAgICAgICAic2N0X0NEMTkiLCJzY3RfTVM0QTEiLCJzY3RfU1RBQjIiLAogICAgICAgICAgICAgICAgICJzY3RfQUxCLTEiLCJzY3RfQ0Q0Iiwic2N0X0NEOEEiLAogICAgICAgICAgICAgICAgICJzY3RfQ0xFQzRHIiwic2N0X0NENUwiLCJzY3RfQzFRQiIsCiAgICAgICAgICAgICAgICAgInNjdF9BQ1RBMiIsInNjdF9WV0YiLCJzY3RfSUdMTDUtMSIsCiAgICAgICAgICAgICAgICAgInNjdF9DRDY4Iiwic2N0X1hDTDE7WENMMiIsInNjdF9MRUYxIiwKICAgICAgICAgICAgICAgICAic2N0X1JTUE8zIiwic2N0X01FQ09NIikKICBnZW5lcyA8LSBjKCJQVFBSQyIsIkNBTENSTCIsIk5LRzciLCJDRDNFIiwiTUFSQ08iLCJMWVoiLAogICAgICAgICAgICAgIkNEMTkiLCJNUzRBMSIsIlNUQUIyIiwiQUxCIiwiQ0Q0IiwiQ0Q4QSIsCiAgICAgICAgICAgICAiQ0xFQzRHIiwiQ0Q1TCIsIkMxUUIiLCJBQ1RBMiIsIlZXRiIsIklHTEw1IiwiQ0Q2OCIsCiAgICAgICAgICAgICAiWENMMTtYQ0wyIiwiTEVGMSIsIlJTUE8zIiwiTUVDT00iKSAKfQpmb3IobnVtIGluIDI6bGVuZ3RoKGdlbmVDb2RlcykpIHsKICBwbG90IDwtIEZlYXR1cmVQbG90KHNvYmosCiAgICAgICAgICAgICAgICAgICAgICBmZWF0dXJlcyA9IGdlbmVDb2Rlc1tudW1dKSArCiAgICBnZ3RpdGxlKGJxdW90ZSh+aXRhbGljKC4oZ2VuZXNbbnVtXSkpKSkKICBwcmludChwbG90KQogIHBkZihwYXN0ZSgiLi9maWd1cmVzLyIsIHRpc3N1ZSwgIi8iLCB0aXNzdWUsICJfIiwgZ2VuZXNbbnVtXSwgIl9VTUFQLnBkZiIsIHNlcCA9ICIiKSkKICBwcmludChwbG90KQogIGRldi5vZmYoKQp9CmBgYAoKIyMjIENlbGwtcG9wdWxhdGlvbi1zcGVjaWZpYyBtYXJrZXJzCgpQbG90IG1hcmtlciBnZW5lcyBmb3Igc3BlY2lmaWMgY2VsbCBwb3B1bGF0aW9ucwoKRW5kb3RoZWxpYWwgY2VsbHM6CgpgYGB7cn0KZ2VuZUNvZGVzIDwtIGMoIlNUQUIyIiwiSVRHQTEiLCJDRDU1IiwiTFlWRTEiLAogICAgICAgICAgICAgICAiQ0QzNCIsIlZXRiIsIklGSVRNMztJRklUTTI7SUZJVE0xIiwiUlNQTzMiLAogICAgICAgICAgICAgICAiTUVDT00iLCJNZWNvbSIsIkNBTENSTCIsIkxPQzExNDA4OTY1NCIsCiAgICAgICAgICAgICAgICJSQU1QMiIsIkJTVDIiLCJDTEVDNEciLCJDVFNWO0NUU0wiLAogICAgICAgICAgICAgICAiU1RBQjEiLCJQTEFDOSIsIlBFQ0FNMSIsIldOVDIiLAogICAgICAgICAgICAgICAiTVlMNiIsIk5SMkYyIikKZ2VuZXMgPC0gYygiU1RBQjIiLCJJVEdBMSIsIkNENTUiLCJMWVZFMSIsCiAgICAgICAgICAgIkNEMzQiLCJWV0YiLCJJRklUTTMiLCJSU1BPMyIsCiAgICAgICAgICAgIk1FQ09NIiwiTWVjb20iLCJDQUxDUkwiLCJMT0MxMTQwODk2NTQiLAogICAgICAgICAgICJSQU1QMiIsIkJTVDIiLCJDTEVDNEciLCJDVFNWO0NUU0wiLAogICAgICAgICAgICJTVEFCMSIsIlBMQUM5IiwiUEVDQU0xIiwiV05UMiIsCiAgICAgICAgICAgIk1ZTDYiLCJOUjJGMiIpIApmb3IobnVtIGluIDE6bGVuZ3RoKGdlbmVDb2RlcykpIHsKICBwbG90IDwtIEZlYXR1cmVQbG90X3NjQ3VzdG9tKHNvYmosCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICBmZWF0dXJlcyA9IGdlbmVDb2Rlc1tudW1dLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgY29sb3JzX3VzZSA9IHZpcmlkaXNfbGlnaHRfaGlnaCkgKwogICAgZ2d0aXRsZShicXVvdGUofml0YWxpYyguKGdlbmVzW251bV0pKSkpCiAgcHJpbnQocGxvdCkKICBwZGYocGFzdGUoIi4vZmlndXJlcy8iLCB0aXNzdWUsICIvIiwgdGlzc3VlLCAiXyIsIGdlbmVzW251bV0sICJfVU1BUC5wZGYiLCBzZXAgPSAiIikpCiAgcHJpbnQocGxvdCkKICBkZXYub2ZmKCkKfQpgYGAKCkx5bXBob2N5dGVzOgoKYGBge3J9CmdlbmVDb2RlcyA8LSBjKCJDWENSMyIsIklMN1ItMSIsCiAgICAgICAgICAgICAgICJTMTAwQTQiLCJHQVRBMyIsIkNEMjQ3IiwiTUFGIiwKICAgICAgICAgICAgICAgIkNDUjQiKQpnZW5lcyA8LSBnZW5lQ29kZXMKZ2VuZUNvZGVzIDwtIGMoIkNEM0QiLCJDRDNFIiwiQ0Q4QSIsIk5LRzciLAogICAgICAgICAgICAgICAiS0xSRDEiLCJLTFJCMSIsIkdaTUEiLCJMRUYxIiwKICAgICAgICAgICAgICAgIklMN1ItMSIsIkVPTUVTIiwiVElHSVQiLCJYQ0wxO1hDTDIiLAogICAgICAgICAgICAgICAiVE9YIiwiQ1RMQTQiLCJDRDQiLCJHWk1LIikKZ2VuZXMgPC0gYygiQ0QzRCIsIkNEM0UiLCJDRDhBIiwiTktHNyIsCiAgICAgICAgICAgIktMUkQxIiwiS0xSQjEiLCJHWk1BIiwiTEVGMSIsCiAgICAgICAgICAgIklMN1IiLCJFT01FUyIsIlRJR0lUIiwiWENMMTtYQ0wyIiwKICAgICAgICAgICAiVE9YIiwiQ1RMQTQiLCJDRDQiLCJHWk1LIikgCmZvcihudW0gaW4gMTpsZW5ndGgoZ2VuZUNvZGVzKSkgewogIHBsb3QgPC0gRmVhdHVyZVBsb3Rfc2NDdXN0b20oc29iaiwKICAgICAgICAgICAgICAgICAgICAgIGZlYXR1cmVzID0gZ2VuZUNvZGVzW251bV0sCiAgICAgICAgICAgICAgICAgICAgICBjb2xvcnNfdXNlID0gdmlyaWRpc19saWdodF9oaWdoKSArCiAgICBnZ3RpdGxlKGJxdW90ZSh+aXRhbGljKC4oZ2VuZXNbbnVtXSkpKSkKICBwcmludChwbG90KQogIHBkZihwYXN0ZSgiLi9maWd1cmVzLyIsIHRpc3N1ZSwgIi8iLCB0aXNzdWUsICJfIiwgZ2VuZXNbbnVtXSwgIl9VTUFQLnBkZiIsIHNlcCA9ICIiKSkKICBwcmludChwbG90KQogIGRldi5vZmYoKQp9CmBgYAoKTHltcGhvY3l0ZXMgd2VpcmQgZ2VuZXM6CgpgYGB7cn0KZ2VuZUNvZGVzIDwtIGMoInNjdF9DQ1I3IikKZ2VuZXMgPC0gZ2VuZUNvZGVzCiNwbG90IDwtIEZlYXR1cmVQbG90X3NjQ3VzdG9tKHNvYmosCiMgICAgICAgICAgICAgICAgICAgICAgZmVhdHVyZXMgPSBnZW5lQ29kZXMsCiMgICAgICAgICAgICAgICAgICAgICAgY29sb3JzX3VzZSA9IHZpcmlkaXNfbGlnaHRfaGlnaCkgKwojICAgIGdndGl0bGUoYnF1b3RlKH5pdGFsaWMoLihnZW5lcykpKSkKcGxvdCA8LSBGZWF0dXJlUGxvdChzb2JqLCBmZWF0dXJlcyA9ICJzY3RfQ0NSNyIpCnByaW50KHBsb3QpCmBgYAoKCk1lc2VuY2h5bWU6CgpgYGB7cn0KZ2VuZUNvZGVzIDwtIGMoIkhISVAiLCJDT0wxQTIiLCJDT0wzQTEiLCJJR0ZCUDciLAogICAgICAgICAgICAgICAiSUdGQlAzIiwiRENOIiwiQ09MMUExIiwiU1BBUkMiLAogICAgICAgICAgICAgICAiUkJQMSIsICJDQUxDUkwiKQpnZW5lcyA8LSBjKCJISElQIiwiQ09MMUEyIiwiQ09MM0ExIiwiSUdGQlA3IiwKICAgICAgICAgICAiSUdGQlAzIiwiRENOIiwiQ09MMUExIiwiU1BBUkMiLAogICAgICAgICAgICJSQlAxIiwiQ0FMQ1JMIikgCmZvcihudW0gaW4gMTpsZW5ndGgoZ2VuZUNvZGVzKSkgewogIHBsb3QgPC0gRmVhdHVyZVBsb3Rfc2NDdXN0b20oc29iaiwKICAgICAgICAgICAgICAgICAgICAgIGZlYXR1cmVzID0gZ2VuZUNvZGVzW251bV0sCiAgICAgICAgICAgICAgICAgICAgICBjb2xvcnNfdXNlID0gdmlyaWRpc19saWdodF9oaWdoKSArCiAgICBnZ3RpdGxlKGJxdW90ZSh+aXRhbGljKC4oZ2VuZXNbbnVtXSkpKSkKICBwcmludChwbG90KQogIHBkZihwYXN0ZSgiLi9maWd1cmVzLyIsIHRpc3N1ZSwgIi8iLCB0aXNzdWUsICJfIiwgZ2VuZXNbbnVtXSwgIl9VTUFQLnBkZiIsIHNlcCA9ICIiKSkKICBwcmludChwbG90KQogIGRldi5vZmYoKQp9CmBgYAoKQ2hvbGFuZ2lvY3l0ZXM6CgpgYGB7cn0KZ2VuZUNvZGVzIDwtIGMoIktSVDE5IiwiQ0ZUUiIsIkVQQ0FNIiwiU0xDNEE0IikKZ2VuZXMgPC0gYygiS1JUMTkiLCJDRlRSIiwiRVBDQU0iLCJTTEM0QTQiKSAKZm9yKG51bSBpbiAxOmxlbmd0aChnZW5lQ29kZXMpKSB7CiAgcGxvdCA8LSBGZWF0dXJlUGxvdF9zY0N1c3RvbShzb2JqLAogICAgICAgICAgICAgICAgICAgICAgZmVhdHVyZXMgPSBnZW5lQ29kZXNbbnVtXSwKICAgICAgICAgICAgICAgICAgICAgIGNvbG9yc191c2UgPSB2aXJpZGlzX2xpZ2h0X2hpZ2gpICsKICAgIGdndGl0bGUoYnF1b3RlKH5pdGFsaWMoLihnZW5lc1tudW1dKSkpKQogIHByaW50KHBsb3QpCiAgcGRmKHBhc3RlKCIuL2ZpZ3VyZXMvIiwgdGlzc3VlLCAiLyIsIHRpc3N1ZSwgIl8iLCBnZW5lc1tudW1dLCAiX1VNQVAucGRmIiwgc2VwID0gIiIpKQogIHByaW50KHBsb3QpCiAgZGV2Lm9mZigpCn0KYGBgCgpMb29rIGF0IHpvbmF0aW9uIGdlbmUgc2lnbmF0dXJlczoKCmBgYHtyfQpjdl9nZW5lcyA8LSBjKCJGRVRVQiIsIkhNR0NTMSIsIkNZUDJFMSIsIkdMVUQxO0dMVUQyIiwiQ1lQMUEyIiwKICAgICAgICAgICAgICAiUkdOIiwiSU5NVCIsIkNPTVQiLCJGRFBTIiwiU1BSLTEiKQpjdl9nZW5lcwpwcF9nZW5lcyA8LSBjKCJTYWEyO1NhYTEtMSIsIkhBTVAiLCJBUE9BMSIsIkFQT0MyIiwiQ1JZTDEiLAogICAgICAgICAgICAgICJBTVkxQTtBTVkxQztBTVkxQjtBTVkyQTtBTVkyQiIsIk1ULUFUUDYiLCJVUk9DMSIsCiAgICAgICAgICAgICAgIkFQT0EyIiwiTVQtQ08zIikKcHBfZ2VuZXMKYGBgCgojIyMgQ2VsbC10eXBlIG1hcmtlcnMgYWxsIGluIG9uZSBwbG90CgpIZXBhdG9jeXRlczoKCmBgYHtyfQpnZW5lQ29kZXMgPC0gYygiQ1lQMkUxIiwiRkVUVUIiLCJITUdDUzEiLCJHTFVEMTtHTFVEMiIsCiAgICAgICAgICAgICAgICJDWVAxQTIiLCJBUE9DMiIsIlNhYTE7IFNhYTIiLCAiSEFNUCIsCiAgICAgICAgICAgICAgICJQQ0sxIikKcGxvdCA8LSBGZWF0dXJlUGxvdF9zY0N1c3RvbShzb2JqLAogICAgICAgICAgICAgICAgICAgICBmZWF0dXJlcyA9IGdlbmVDb2RlcywKICAgICAgICAgICAgICAgICAgICAgY29sb3JzX3VzZSA9IHZpcmlkaXNfbGlnaHRfaGlnaCwKICAgICAgICAgICAgICAgICAgICAgbnVtX2NvbHVtbnMgPSA0KQpwbG90CnBkZihwYXN0ZSgiLi9maWd1cmVzLyIsIHRpc3N1ZSwgIi8iLCB0aXNzdWUsICJfaGVwYXRvY3l0ZU1hcmtlcnNfVU1BUC5wZGYiLCBzZXAgPSAiIiksCiAgICBoZWlnaHQgPSBjZWlsaW5nKGxlbmd0aChnZW5lQ29kZXMpLzQpKjYsIHdpZHRoID0gMjQsIHdpZHRoID0gMjQpCnBsb3QKZGV2Lm9mZigpCmBgYAoKRW5kb3RoZWxpYWwgY2VsbHM6CgpgYGB7cn0KIyBPbmUgc2V0IG9mIGZlYXR1cmVzCmdlbmVDb2RlcyA8LSBjKCJTVEFCMiIsIklUR0ExIiwiTFlWRTEiLCJJRklUTTM7SUZJVE0yO0lGSVRNMSIsCiAgICAgICAgICAgICAgICJDQUxDUkwiLCJSQU1QMiIsIkJTVDIiLCAiQ0xFQzRHIiwKICAgICAgICAgICAgICAgIkNUU1Y7Q1RTTCIsICJTVEFCMSIsICJWV0YiLCAiQ0QzNCIsCiAgICAgICAgICAgICAgICJNRUNPTSIsIlBMQUM5IiwiUlNQTzMiLCJXTlQyIikKIyBTZWNvbmQgc2V0IG9mIGZlYXR1cmVzCmdlbmVDb2RlcyA8LSBjKCJDTEVDNEciLCJTVEFCMSIsIlNUQUIyIiwiTFlWRTEiLAogICAgICAgICAgICAgICAiQ0QzNCIsIk1FQ09NIiwiVldGIiwiUlNQTzMiLAogICAgICAgICAgICAgICAiV05UMiIsIkFDS1IxIikKcGxvdCA8LSBGZWF0dXJlUGxvdF9zY0N1c3RvbShzb2JqLAogICAgICAgICAgICAgICAgICAgICBmZWF0dXJlcyA9IGdlbmVDb2RlcywKICAgICAgICAgICAgICAgICAgICAgY29sb3JzX3VzZSA9IHZpcmlkaXNfbGlnaHRfaGlnaCwKICAgICAgICAgICAgICAgICAgICAgbnVtX2NvbHVtbnMgPSA0KQpwbG90CnBkZihwYXN0ZSgiLi9maWd1cmVzLyIsIHRpc3N1ZSwgIi8iLCB0aXNzdWUsICJfZW5kb3RoZWxpYWxNYXJrZXJzX1VNQVAucGRmIiwgc2VwID0gIiIpLAogICAgaGVpZ2h0ID0gY2VpbGluZyhsZW5ndGgoZ2VuZUNvZGVzKS80KSo2LCB3aWR0aCA9IDI0LCB3aWR0aCA9IDI0KQpwbG90CmRldi5vZmYoKQpgYGAKClQgbHltcGhvY3l0ZXMKCmBgYHtyfQpnZW5lQ29kZXMgPC0gYygiQ0QzRCIsIkNEM0UiLCJDRDQiLCJDRDhBIiwKICAgICAgICAgICAgICAgIkNUTEE0IiwiSUw3Ui0xIiwiTEVGMSIsIkVPTUVTIiwKICAgICAgICAgICAgICAgIlRJR0lUIiwiS0xSQjEiLCJLTFJEMSIsIk5LRzciLAogICAgICAgICAgICAgICAiWENMMTtYQ0wyIiwiVE9YIiwiR1pNQSIsIkdaTUsiKQpnZW5lQ29kZXMgPC0gYygiQ0QzRCIsIkNEM0UiLCJHSU1BUDEiLCJHSU1BUDQiLAogICAgICAgICAgICAgICAiR0lNQVA2IiwiR0lNQVA3IiwiQ0Q0IiwiQ0QyOCIsCiAgICAgICAgICAgICAgICJGT1hQMSIsIkNUTEE0IiwiQ0QzOCIsIlhDTDE7WENMMiIsCiAgICAgICAgICAgICAgICJJTDJSQSIsIk5LRzciLCJDRDhBIiwiQ0NMNSIsCiAgICAgICAgICAgICAgICJDRDY5IiwiQ0QyIiwiR1pNSyIsIkJSQ0EyIiwgCiAgICAgICAgICAgICAgICJUT1AyQSIsICJTVE1OMSIpCnBsb3QgPC0gRmVhdHVyZVBsb3Rfc2NDdXN0b20oc29iaiwKICAgICAgICAgICAgICAgICAgICAgZmVhdHVyZXMgPSBnZW5lQ29kZXMsCiAgICAgICAgICAgICAgICAgICAgIGNvbG9yc191c2UgPSB2aXJpZGlzX2xpZ2h0X2hpZ2gsCiAgICAgICAgICAgICAgICAgICAgIG51bV9jb2x1bW5zID0gNCkKcGxvdApwZGYocGFzdGUoIi4vZmlndXJlcy8iLCB0aXNzdWUsICIvIiwgdGlzc3VlLCAiX2x5bXBob2N5dGVNYXJrZXJzX1VNQVAucGRmIiwgc2VwID0gIiIpLAogICAgaGVpZ2h0ID0gY2VpbGluZyhsZW5ndGgoZ2VuZUNvZGVzKS80KSo2LCB3aWR0aCA9IDI0KQpwbG90CmRldi5vZmYoKQpgYGAKClQgbHltcGhvY3l0ZXMgZm9yIFBCTUNzCgpgYGB7cn0KZ2VuZUNvZGVzIDwtIGMoIkNEM0QiLCJDRDNFIiwiWENMMTtYQ0wyIiwiTktHNyIsCiAgICAgICAgICAgICAgICJMRUYxIiwiSUw3Ui0xIiwiR0FUQTMiLCJNQUYiLAogICAgICAgICAgICAgICAiQ0NSNCIsIlMxMDBBNCIsIkNYQ1IzIiwiSUwyUkEiLAogICAgICAgICAgICAgICAiQ0QxNjAiLCJDRDI0NyIsIkJSQ0EyIiwiU1RNTjEiKQojIENDUjcgYW5kIFRDRjcgYXJlbid0IHNob3dpbmcgdXA/Pz8KcGxvdCA8LSBGZWF0dXJlUGxvdF9zY0N1c3RvbShzb2JqLAogICAgICAgICAgICAgICAgICAgICBmZWF0dXJlcyA9IGdlbmVDb2RlcywKICAgICAgICAgICAgICAgICAgICAgY29sb3JzX3VzZSA9IHZpcmlkaXNfbGlnaHRfaGlnaCwKICAgICAgICAgICAgICAgICAgICAgbnVtX2NvbHVtbnMgPSA0LAogICAgICAgICAgICAgICAgICAgICBzbG90ID0gImRhdGEiKQpwbG90CnBkZihwYXN0ZSgiLi9maWd1cmVzLyIsIHRpc3N1ZSwgIi8iLCB0aXNzdWUsICJfbHltcGhvY3l0ZU1hcmtlcnNfVU1BUC5wZGYiLCBzZXAgPSAiIiksCiAgICBoZWlnaHQgPSBjZWlsaW5nKGxlbmd0aChnZW5lQ29kZXMpLzQpKjYsIHdpZHRoID0gMjQpCnBsb3QKZGV2Lm9mZigpCmBgYAoKCk1lc2VuY2h5bWU6CgpgYGB7cn0KZ2VuZUNvZGVzIDwtIGMoIkRDTiIsIkNPTDFBMiIsIkNPTDNBMSIsIkhISVAiLAogICAgICAgICAgICAgICAiQ0FMQ1JMIiwiU1BBUkMiLCJSQlAxIiwiQ09MMUExIiwKICAgICAgICAgICAgICAgIklHRkJQMyIsICJJR0ZCUDciKQpnZW5lQ29kZXMgPC0gYygiQ09MMUEyIiwgIkNPTDNBMSIsICJJR0ZCUDciLCAiSUdGQlAzIiwgCiAgICAgICAgICAgICAgICJEQ04iLCAiQ0FMQ1JMIiwgIkNPTDFBMSIsICJTUEFSQyIsIAogICAgICAgICAgICAgICAiUkJQMSIsICJISElQIiwiUkJQMSIsIkxSQVQiLAogICAgICAgICAgICAgICAiUERFM0IiLCJIR0YiLCJDTk4yIiwiQUNUQTIiKQpwbG90IDwtIEZlYXR1cmVQbG90X3NjQ3VzdG9tKHNvYmosCiAgICAgICAgICAgICAgICAgICAgIGZlYXR1cmVzID0gZ2VuZUNvZGVzLAogICAgICAgICAgICAgICAgICAgICBjb2xvcnNfdXNlID0gdmlyaWRpc19saWdodF9oaWdoLAogICAgICAgICAgICAgICAgICAgICBudW1fY29sdW1ucyA9IDQpCnBsb3QKcGRmKHBhc3RlKCIuL2ZpZ3VyZXMvIiwgdGlzc3VlLCAiLyIsIHRpc3N1ZSwgIl9tZXNlbmNoeW1hbE1hcmtlcnNfVU1BUC5wZGYiLCBzZXAgPSAiIiksCiAgICBoZWlnaHQgPSBjZWlsaW5nKGxlbmd0aChnZW5lQ29kZXMpLzQpKjYsIHdpZHRoID0gMjQpCnBsb3QKZGV2Lm9mZigpCmBgYAoKQ2hvbGFuZ2lvY3l0ZXM6CgpgYGB7cn0KZ2VuZUNvZGVzIDwtIGMoIktSVDE5IiwiU0xDNEE0IiwiQ0ZUUiIsIkVQQ0FNIikKZ2VuZUNvZGVzIDwtIGMoIkFOWEEyIiwgIkNTVDMiLCAiQklSQzMiLCAiVEVTQyIsCiAgICAgICAgICAgICAgICJLUlQxOSIsICJTT1g5IiwgIkVQQ0FNIiwiU0xDNEE0IiwKICAgICAgICAgICAgICAgIkNGVFIiKQpwbG90IDwtIEZlYXR1cmVQbG90X3NjQ3VzdG9tKHNvYmosCiAgICAgICAgICAgICAgICAgICAgIGZlYXR1cmVzID0gZ2VuZUNvZGVzLAogICAgICAgICAgICAgICAgICAgICBjb2xvcnNfdXNlID0gdmlyaWRpc19saWdodF9oaWdoLAogICAgICAgICAgICAgICAgICAgICBudW1fY29sdW1ucyA9IDQpCnBsb3QKcGRmKHBhc3RlKCIuL2ZpZ3VyZXMvIiwgdGlzc3VlLCAiLyIsIHRpc3N1ZSwgIl9jaG9sYW5naW9jeXRlTWFya2Vyc19VTUFQLnBkZiIsIHNlcCA9ICIiKSwKICAgIGhlaWdodCA9IGNlaWxpbmcobGVuZ3RoKGdlbmVDb2RlcykvNCkqNiwgd2lkdGggPSAyNCkKcGxvdApkZXYub2ZmKCkKYGBgCgpNeWVsb2lkIGNlbGxzCgpgYGB7cn0KZ2VuZUNvZGVzIDwtIGMoIlRZUk9CUCIsIkNENzQiLCJDVFNTIiwiVlNJRzQiLAogICAgICAgICAgICAgICAiTUFSQ08iLCAiQ0Q1TCIsIkhNT1gxIiwiVkNBTTEiLAogICAgICAgICAgICAgICAiSVRHQU0iLCJGQ0dSM0E7RkNHUjNCIiwiRmNncjM7RmNncjJiIiwiQ1RTQSIsCiAgICAgICAgICAgICAgICJJTDE3UkEtMSIsIlNJUlBBO1NJUlBCMTtTSVJQRyIsIkJUTEEiLCAiQ0FETTEiLAogICAgICAgICAgICAgICAiSUQyIiwiTFk2RSIsIkZMVDMiLCAiQ0QxNCIsCiAgICAgICAgICAgICAgICJOT1RDSDIiLCJMWVotMSIsIlMxMDBBOCIpCnBsb3QgPC0gRmVhdHVyZVBsb3Rfc2NDdXN0b20oc29iaiwKICAgICAgICAgICAgICAgICAgICAgZmVhdHVyZXMgPSBnZW5lQ29kZXMsCiAgICAgICAgICAgICAgICAgICAgIGNvbG9yc191c2UgPSB2aXJpZGlzX2xpZ2h0X2hpZ2gsCiAgICAgICAgICAgICAgICAgICAgIG51bV9jb2x1bW5zID0gNCkKcGxvdApwZGYocGFzdGUoIi4vZmlndXJlcy8iLCB0aXNzdWUsICIvIiwgdGlzc3VlLCAiX215ZWxvaWRNYXJrZXJzX1VNQVAucGRmIiwgc2VwID0gIiIpLAogICAgaGVpZ2h0ID0gY2VpbGluZyhsZW5ndGgoZ2VuZUNvZGVzKS80KSo2LCB3aWR0aCA9IDI0KQpwbG90CmRldi5vZmYoKQpgYGAKCgojIyBDb3JyZWxhdGlvbiBwbG90cwoKR2VuZXJhdGUgaGVhdG1hcHMgY29tcGFyaW5nIHdvb2RjaHVjayBjbHVzdGVycyB3aXRoIHZhcmlvdXMgZGF0YXNldHMKCiMjIyBTZXR1cCBmb3IgYWxsIGNvcnJlbGF0aW9ucwoKRmlyc3QsIHNldCB1cCB3aGF0ZXZlciB3b29kY2h1Y2sgZGF0YXNldCBJIGFtIHdvcmtpbmcgd2l0aCBieSByZWFkaW5nIGluIHRoZSBvcnRob2xvZyB0YWJsZSwgY2hvb3Npbmcgd2hldGhlciB0aGUgb3J0aG9sb2dzIHRvIHVzZSBhcmUgaHVtYW4sIG1vdXNlLCBvciB3b29kY2h1Y2ssIGFuZCBjYWxjdWxhdGluZyB0aGUgYXZlcmFnZSBleHByZXNzaW9uIGZvciBlYWNoIGNsdXN0ZXIuIFRoZSBvdXRwdXQgb2YgdGhpcyBzZWN0aW9uIGlzIHRoZSBzY2FsZWQgY2x1c3RlciBnZW5lLWV4cHJlc3Npb24gbWF0cml4CgpgYGB7cn0KI2dyb3VwcyA8LSByZXMKZ3JvdXBzIDwtICJnZW5lcmFsX2NlbGxfbGFiZWxzIgojIFJlYWQgaW4gb3J0aG9sb2cgdGFibGUKZ2VuZU5hbWVUYWJsZSA8LSByZWFkLnRhYmxlKCJ+L0Ryb3Bib3gvWm9lL3NjZl92ZXJzaW9uL21ha2VfZ3RmL29ydGhvZmluZGVyX3NjMi9ob21vbG9nZW5lL2NvbGxlY3RlZE9ydGhvZmluZGVyUGFpcmluZ3MudHN2IiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgIHNlcCA9ICJcdCIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICBoZWFkZXIgPSBUUlVFKQp3b29kY2h1Y2tDbHVzdGVyQXZlcmFnZXMgPC0gQXZlcmFnZUV4cHJlc3Npb24oc29iaiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIGFzc2F5cyA9ICJTQ1QiLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgc2xvdCA9ICJzY2FsZS5kYXRhIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIGdyb3VwLmJ5ID0gZ3JvdXBzKQojIFNjYWxlIGRhdGEKd29vZGNodWNrQ2x1c3RlckF2ZXJhZ2VzJFNDVCA8LSBuYS5vbWl0KHQoc2NhbGUodChhcy5tYXRyaXgod29vZGNodWNrQ2x1c3RlckF2ZXJhZ2VzJFNDVCkpKSkpCiMgR3JhYiBnZW5lIG5hbWVzIGZyb20gU2V1cmF0IG9iamVjdAp1bmlxdWVIaWVyIDwtIHJvdy5uYW1lcyh3b29kY2h1Y2tDbHVzdGVyQXZlcmFnZXMkU0NUKQp1bmlxdWVIaWVyIDwtIGFzLmRhdGEuZnJhbWUodW5pcXVlSGllcikKIyBCaW5kIHdpdGggZ2VuZU5hbWVUYWJsZSB0byBnZXQgY29ycmVjdCBvcmRlciAobm90aWNlIHVuaXF1ZUhpZXIgaXMgb24gbGVmdCkKbmV3TmFtZXMgPC0gZHBseXI6OmxlZnRfam9pbih1bmlxdWVIaWVyLCBnZW5lTmFtZVRhYmxlLCBieSA9ICJ1bmlxdWVIaWVyIikKIyBHZXQgb3J0aG9sb2dzIGZyb20gZWl0aGVyIG1vdXNlIG9yIGh1bWFuCnNwZWNpZXMgPC0gImh1bWFuIgppZiAoc3BlY2llcyA9PSAiaHVtYW4iKSB7CiAgIyBJZiBodW1hbiBvbmUtdG8tb25lIG9ydGhvIGhhcyBOQSwgcmVwbGFjZSB3aXRoIG1pa2Fkb19maW5hbF9zYzJfc3RyaW5nZW50X25vTWl0b19wcm90ZWluIGNvbHVtbgogICMgVGhpcyBpcyB0byBhdm9pZCBhbmQgcG90ZW50aWFsIG1pc3Rha2VzIGluIHJlY29nbml6aW5nIHRoaW5ncyBpdCBzaG91bGRuJ3QgYmUgcmVjb2duaXppbmcKICBuZXdOYW1lcyRzcGVjaWVzT25lVG9PbmUgPC0gaWZlbHNlKGlzLm5hKG5ld05hbWVzJGh1bWFuT25lVG9PbmUpLCBuZXdOYW1lcyR1bmlxdWVIaWVyLCBuZXdOYW1lcyRodW1hbk9uZVRvT25lKQp9IGVsc2UgaWYgKHNwZWNpZXMgPT0gIm1vdXNlIikgewogIG5ld05hbWVzJHNwZWNpZXNPbmVUb09uZSA8LSBpZmVsc2UoaXMubmEobmV3TmFtZXMkbW91c2VPbmVUb09uZSksIG5ld05hbWVzJHVuaXF1ZUhpZXIsIG5ld05hbWVzJG1vdXNlT25lVG9PbmUpCn0gZWxzZSBpZiAoc3BlY2llcyA9PSAid29vZGNodWNrIikgewogIG5ld05hbWVzJHNwZWNpZXNPbmVUb09uZSA8LSBuZXdOYW1lcyR1bmlxdWVIaWVyCn0KIyBHcmFiIGRhdGFmcmFtZQp3b29kY2h1Y2tDbHVzdGVyQXZlcmFnZXMgPC0gd29vZGNodWNrQ2x1c3RlckF2ZXJhZ2VzJFNDVAojIFJlcGxhY2UgbmFtZXMgd2l0aCBvbmUtdG8tb25lIG9ydGhvbG9ndWUgb2YgcGFydGljdWxhciBzcGVjaWVzCnJvdy5uYW1lcyh3b29kY2h1Y2tDbHVzdGVyQXZlcmFnZXMpIDwtIG5ld05hbWVzJHNwZWNpZXNPbmVUb09uZQojIE1ha2Ugc3VyZSBmb3JtYXR0ZWQgY29ycmVjdGx5Cndvb2RjaHVja0NsdXN0ZXJBdmVyYWdlcyA8LSBhcy5kYXRhLmZyYW1lKHdvb2RjaHVja0NsdXN0ZXJBdmVyYWdlcykKIyBPcmRlciBieSBnZW5lIG5hbWUKd29vZGNodWNrQ2x1c3RlckF2ZXJhZ2VzIDwtIHdvb2RjaHVja0NsdXN0ZXJBdmVyYWdlc1tvcmRlcihyb3cubmFtZXMod29vZGNodWNrQ2x1c3RlckF2ZXJhZ2VzKSksXQojIFNhbml0eSBjaGVjawpoZWFkKHdvb2RjaHVja0NsdXN0ZXJBdmVyYWdlcykKYGBgCgojIyMgRGF0YXNldC1zcGVjaWZpYyBzZXR1cAoKQ29ycmVsYXRpb24gb2Ygd29vZGNodWNrIFBCTUNzIHdpdGggaHVtYW4gNjhrIFBCTUMgZGF0YXNldCBmcm9tIDEwWCBHZW5vbWljcwoKYGBge3J9CiMgTm93IG5lZWQgdG8gcmVhZCBpbiA2OEsgUEJNQyBkYXRhCmh1bWFuUEJNQyA8LSByZWFkLmNzdigifi9Ecm9wYm94L1pvZS9zY2ZfdmVyc2lvbi9hbmFseXNpcy9jb3JyZWxhdGlvblRlc3RzLzY4S19wYm1jX2RhdGEvNjhLX2VucmljaGVkR2VuZXMuY3N2IiwKICAgICAgICAgICAgICAgICAgICAgIGhlYWRlciA9IEZBTFNFKQojIEdldCByaWQgb2YgdG9wIHJvdwpodW1hblBCTUMgPC0gaHVtYW5QQk1DWy0xLF0KIyBTZXBhcmF0ZSBhbGwgY2VsbCBhbmQgbXlsb2lkIGNlbGwgZGF0YQphbGxDZWxscyA8LSBzZWxlY3QoaHVtYW5QQk1DLCBWMSwgVjIsIFYzKQpteWVsb2lkIDwtIHNlbGVjdChodW1hblBCTUMsIFY1LCBWNiwgVjcpCiMgUmVuYW1lIGFuZCBnZXQgcmlkIG9mIGZpcnN0IHJvdwpjb2xuYW1lcyhhbGxDZWxscykgPC0gYygiQ2x1c3RlciIsICJHZW5lIiwgIkVucmljaG1lbnQiKQpjb2xuYW1lcyhteWVsb2lkKSA8LSBjKCJDbHVzdGVyIiwgIkdlbmUiLCAiRW5yaWNobWVudCIpCiMgR2V0IHJpZCBvZiB0b3Agcm93CmFsbENlbGxzIDwtIGFsbENlbGxzWy0xLF0KbXllbG9pZCA8LSBteWVsb2lkWy0xLF0KIyBJIGNhbiB0aGVuIGZpbHRlciB0aGUgcm93cyBhbmQgYmluZCB0aGUgZGF0YSBmcmFtZXMgYmFjayB0b2dldGhlciBieSBnZW5lIG5hbWUKZm9yIChqIGluIDE6MTApIHsKICBkYXQgPC0gZHBseXI6OmZpbHRlcihhbGxDZWxscywgQ2x1c3RlciA9PSBqKQogICNkYXQgPC0gZ2V0KHBhc3RlKCJhbGxDZWxscyIsIGosIHNlcCA9ICIiKSkKICAjIE11bHRpcGx5IGVucmljaG1lbnQgdmFsdWVzIGJ5IC0xIGJlY2F1c2UgdGhlIHNpZ25zIGFyZSBiYWNrd2FyZHM/Pz8KICBkYXQkRW5yaWNobWVudCA8LSBhcy5udW1lcmljKGRhdCRFbnJpY2htZW50KSAqIC0xCiAgZGF0IDwtIGRwbHlyOjpzZWxlY3QoZGF0LCBHZW5lLCBFbnJpY2htZW50KQogIGNvbG5hbWVzKGRhdCkgPC0gYygiR2VuZSIsIHBhc3RlKCJFbnJpY2htZW50Iiwgaiwgc2VwID0gIiIpKQogIGlmIChqID09IDEpIHsKICAgIGFsbENlbGxzTWF0cml4IDwtIGRhdAogIH0KICBlbHNlIHsKICAgIGFsbENlbGxzTWF0cml4IDwtIGRwbHlyOjpmdWxsX2pvaW4oYWxsQ2VsbHNNYXRyaXgsIGRhdCwgYnkgPSAiR2VuZSIpCiAgfQp9CiMgTWFrZSByb3cgbmFtZXMgZ2VuZSBuYW1lcwpyb3duYW1lcyhhbGxDZWxsc01hdHJpeCkgPC0gYWxsQ2VsbHNNYXRyaXgkR2VuZQphbGxDZWxsc01hdHJpeCA8LSBkcGx5cjo6c2VsZWN0KGFsbENlbGxzTWF0cml4LCAtR2VuZSkKIyBNYWtlIGNvbHVtbiBuYW1lcyBjZWxsIHR5cGVzCmNvbG5hbWVzKGFsbENlbGxzTWF0cml4KSA8LSBjKCJBY3RpdmF0ZWQgQ0Q4KyIsICJOYWl2ZSBDRDgrIiwgIk1lbW9yeSBhbmQgUmVnIFQiLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAiTmFpdmUgQ0Q0KyIsICJOSyIsICJDRDgrIiwgIkIiLCAiTWVnYWthcnlvY3l0ZXMiLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAiTW9ub2N5dGVzIGFuZCBEZW5kcml0aWMiLCAiQiwgRGVuZHJpdGljLCBUIikKIyBTY2FsZSBhY3Jvc3MgY29sdW1ucwphbGxDZWxsc01hdHJpeCA8LSB0KHNjYWxlKHQoYWxsQ2VsbHNNYXRyaXgpKSkKIyBPcmRlciBnZW5lcyBhbHBoYWJldGljYWxseSBieSBnZW5lIG5hbWUKYWxsQ2VsbHNNYXRyaXggPC0gYWxsQ2VsbHNNYXRyaXhbb3JkZXIocm93Lm5hbWVzKGFsbENlbGxzTWF0cml4KSksXQojIFNhbml0eSBjaGVjawpoZWFkKGFsbENlbGxzTWF0cml4KQpzcGVjaWVzRGF0YSA8LSAiNjhLUEJNQyIKYGBgCgpDb3JyZWxhdGlvbiBvZiB3b29kY2h1Y2sgbGl2ZXIgd2l0aCBodW1hbiBsaXZlciBkYXRhc2V0IGZyb20gTWFjUGFybGFuZCAqZXQgYWwuKiAoMjAxOCkKCmBgYHtyfQojIEZpbmQgY2x1c3RlciBhdmVyYWdlcyBvZiBodW1hbiBsaXZlciBkYXRhCmxvYWQoIn4vRHJvcGJveC9ab2Uvc2NmX3ZlcnNpb24vYW5hbHlzaXMvY29ycmVsYXRpb25UZXN0cy9IdW1hbkxpdmVyLlJEYXRhIikKIyBSdW4gU0NUcmFuc2Zvcm0KSHVtYW5MaXZlclNldXJhdCA8LSBVcGRhdGVTZXVyYXRPYmplY3QoSHVtYW5MaXZlclNldXJhdCkKSHVtYW5MaXZlclNldXJhdCA8LSBTQ1RyYW5zZm9ybShIdW1hbkxpdmVyU2V1cmF0KQpodW1hbkNsdXN0ZXJBdmVyYWdlcyA8LSBBdmVyYWdlRXhwcmVzc2lvbihIdW1hbkxpdmVyU2V1cmF0LAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICBhc3NheXMgPSAiU0NUIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgc2xvdCA9ICJzY2FsZS5kYXRhIikKIyBSZXBsYWNlIGNsdXN0ZXIgbnVtYmVycyB3aXRoIG5hbWVzCmNvbG5hbWVzKGh1bWFuQ2x1c3RlckF2ZXJhZ2VzJFNDVCkgPC0gYygiSGVwIDEiLCAiQWxwaGEtYmV0YSBUIGNlbGxzIiwgIkhlcCAyIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICJJbmZsYW1tYXRvcnkgbWFjcyIsICJIZXAgMyIsICJIZXAgNCIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAiUGxhc21hIGNlbGxzIiwgIk5LLWxpa2UgY2VsbHMiLCAiR2FtbWEtZGVsdGEgVCBjZWxscyIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAiTm9uLWluZmxhbW1hdG9yeSBtYWNzIiwgIlBlcmlwb3J0YWwgTFNFQ3MiLCAiQ2VudHJhbCB2ZW5vdXMgTFNFQ3MiLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIlBvcnRhbCBlbmRvdGhlbGlhbCBjZWxscyIsICJIZXAgNSIsICJIZXAgNiIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAiTWF0dXJlIEIgY2VsbHMiLCAiQ2hvbGFuZ2lvY3l0ZXMiLCAiR2FtbWEtZGVsdGEgVCBjZWxscyAyIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICJFcnl0aHJvaWQgY2VsbHMiLCAiSGVwYXRpYyBzdGVsbGF0ZSBjZWxscyIpCiMgSWYgb25seSBsb29raW5nIGF0IHNwZWNpZmljIGNsdXN0ZXJzCiNodW1hbkNsdXN0ZXJBdmVyYWdlcyRTQ1QgPC0gaHVtYW5DbHVzdGVyQXZlcmFnZXMkU0NUWyxjKCIzIiwiMSIsIjE1IiwiNiIsIjE0IiwiNSIpXQojIE90aGVyd2lzZSBnbyBzdHJhaWdodCB0byBoZXJlOgpodW1hbkNsdXN0ZXJBdmVyYWdlcyRTQ1QgPC0gbmEub21pdCh0KHNjYWxlKHQoYXMubWF0cml4KGh1bWFuQ2x1c3RlckF2ZXJhZ2VzJFNDVCkpKSkpCiMgR3JhYiBnZW5lIG5hbWVzCmh1bWFuR2VuZXMgPC0gcm93Lm5hbWVzKGh1bWFuQ2x1c3RlckF2ZXJhZ2VzJFNDVCkKIyBOb3cgdHVybiBpbnRvIGxhcmdlIGRhdGFmcmFtZQphbGxDZWxsc01hdHJpeCA8LSBhcy5kYXRhLmZyYW1lKGh1bWFuQ2x1c3RlckF2ZXJhZ2VzJFNDVCkKIyBPcmRlciBieSByb3cgbmFtZQphbGxDZWxsc01hdHJpeCA8LSBhbGxDZWxsc01hdHJpeFtvcmRlcihyb3cubmFtZXMoYWxsQ2VsbHNNYXRyaXgpKSxdCiMgU2FuaXR5IGNoZWNrCmhlYWQoYWxsQ2VsbHNNYXRyaXgpCnNwZWNpZXNEYXRhIDwtICJtYWNwYXJsYW5kIgpgYGAKCkNvcnJlbGF0aW9uIG9mIHdvb2RjaHVjayBsaXZlciB3aXRoIGh1bWFuIGxpdmVyIGRhdGFzZXQgZnJvbSBBaXphcmFuaSAqZXQgYWwuKgoKYGBge3J9CiMgUmVhZCBpbiBBaXphcmFuaSBkYXRhc2V0CmFpemFyYW5pIDwtIHJlYWRSRFMoIn4vRHJvcGJveC9ab2Uvc2NmX3ZlcnNpb24vYW5hbHlzaXMvY29ycmVsYXRpb25UZXN0cy9HU0UxMjQzOTVfTm9ybWFsaHVtYW5saXZlcmRhdGEuUkRhdGEiKQojIFJlYWQgaW4gY2x1c3RlcnMgYW5kIGxhYmVsIGNlbGxzCmFpemFyYW5pQ2x1c3RlcnMgPC0gcmVhZC50YWJsZSgifi9Ecm9wYm94L1pvZS9zY2ZfdmVyc2lvbi9hbmFseXNpcy9jb3JyZWxhdGlvblRlc3RzL0dTRTEyNDM5NV9jbHVzdGVycGFydGl0aW9uLnR4dCIpCiMgT25seSBrZWVwIGNlbGxzIGluIHRoZSBjbHVzdGVyIG9iamVjdAphaXphcmFuaSA8LSBhaXphcmFuaVssaW50ZXJzZWN0KGNvbG5hbWVzKGFpemFyYW5pKSxyb3cubmFtZXMoYWl6YXJhbmlDbHVzdGVycykpXQojIENyZWF0ZSBTZXVyYXQgb2JqZWN0CmFpemFyYW5pIDwtIENyZWF0ZVNldXJhdE9iamVjdChjb3VudHMgPSBhaXphcmFuaSkKIyBSdW4gU0NUcmFuc2Zvcm0KYWl6YXJhbmkgPC0gU0NUcmFuc2Zvcm0oYWl6YXJhbmkpCiMgQWRkIGNsdXN0ZXIgSURzCklkZW50cyhhaXphcmFuaSkgPC0gYWl6YXJhbmlDbHVzdGVycyRzY3QuY3BhcnQKIyBHZXQgY2x1c3RlciBhdmVyYWdlcwphaXphcmFuaUF2ZXJhZ2VzIDwtIEF2ZXJhZ2VFeHByZXNzaW9uKGFpemFyYW5pLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIGFzc2F5cyA9ICJTQ1QiLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIHNsb3QgPSAic2NhbGUuZGF0YSIpCmFpemFyYW5pQXZlcmFnZXMkU0NUIDwtIG5hLm9taXQodChzY2FsZSh0KGFzLm1hdHJpeChhaXphcmFuaUF2ZXJhZ2VzJFNDVCkpKSkpCiMgR3JhYiBnZW5lIG5hbWVzCmFpemFyYW5pR2VuZXMgPC0gcm93Lm5hbWVzKGFpemFyYW5pQXZlcmFnZXMkU0NUKQojIE5vdyB0dXJuIGludG8gbGFyZ2UgZGF0YWZyYW1lCmFsbENlbGxzTWF0cml4IDwtIGFzLmRhdGEuZnJhbWUoYWl6YXJhbmlBdmVyYWdlcyRTQ1QpCiMgT3JkZXIgYnkgY29sdW1uCmFsbENlbGxzTWF0cml4IDwtIGFsbENlbGxzTWF0cml4Wyxhcy5jaGFyYWN0ZXIoc29ydChhcy5udW1lcmljKGNvbG5hbWVzKGFsbENlbGxzTWF0cml4KSkpKV0KIyBSZW5hbWUgY29sdW1ucyB0byBiZSBtb3JlIG1lYW5pbmdmdWwgKG5vdCB0b3RhbGx5IGNvbmZpZGVudCBJIGdvdCBhbGwgY29ycmVjdCkKY29sbmFtZXMoYWxsQ2VsbHNNYXRyaXgpIDwtIGMoJ05LLCBOS1QgYW5kIFQgY2VsbHMgKDEpJywKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgJ0t1cGZmZXIgY2VsbHMgKDIpJywKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgJ05LLCBOS1QgYW5kIFQgY2VsbHMgKDMpJywKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgJ0VQQ0FNKyBjZWxscyBhbmQgY2hvbGFuZ2lvY3l0ZXMgKDQpJywKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgJ05LLCBOS1QgYW5kIFQgY2VsbHMgKDUpJywKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgJ0t1cGZmZXIgY2VsbHMgKDYpJywKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgJ0VQQ0FNKyBjZWxscyBhbmQgY2hvbGFuZ2lvY3l0ZXMgKDcpJywKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgJ0IgY2VsbHMgKDgpJywKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgJ0xpdmVyIHNpbnVzb2lkYWwgZW5kb3RoZWxpYWwgY2VsbHMgKDkpJywKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgJ01hY3JvdmFzY3VsYXIgZW5kb3RoZWxpYWwgY2VsbHMgKDEwKScsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICdIZXBhdG9jeXRlcyAoMTEpJywKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgJ05LLCBOS1QgYW5kIFQgY2VsbHMgKDEyKScsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICdMaXZlciBzaW51c29pZGFsIGVuZG90aGVsaWFsIGNlbGxzICgxMyknLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAnSGVwYXRvY3l0ZXMgKDE0KScsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICdPdGhlciBlbmRvdGhlbGlhbCBjZWxscyAoMTUpJywKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgJ090aGVyICgxNiknLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAnSGVwYXRvY3l0ZXMgKDE3KScsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICdOSywgTktUIGFuZCBUIGNlbGxzICgxOCknLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAnTkssIE5LVCBhbmQgVCBjZWxscyAoMTkpJywKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgJ0xpdmVyIHNpbnVzb2lkYWwgZW5kb3RoZWxpYWwgY2VsbHMgKDIwKScsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICdNYWNyb3Zhc2N1bGFyIGVuZG90aGVsaWFsIGNlbGxzICgyMSknLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAnQiBjZWxscyAoMjIpJywKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgJ0t1cGZmZXIgY2VsbHMgKDIzKScsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICdFUENBTSsgY2VsbHMgYW5kIGNob2xhbmdpb2N5dGVzICgyNCknLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAnS3VwZmZlciBjZWxscyAoMjUpJywKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgJ090aGVyIGVuZG90aGVsaWFsIGNlbGxzICgyNiknLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAnT3RoZXIgKDI3KScsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICdOSywgTktUIGFuZCBUIGNlbGxzICgyOCknLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAnTWFjcm92YXNjdWxhciBlbmRvdGhlbGlhbCBjZWxscyAoMjkpJywKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgJ0hlcGF0b2N5dGVzICgzMCknLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAnS3VwZmZlciBjZWxscyAoMzEpJywKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgJ01hY3JvdmFzY3VsYXIgZW5kb3RoZWxpYWwgY2VsbHMgKDMyKScsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICdTdGVsbGF0ZSBjZWxscyBhbmQgbXlvZmlicm9ibGFzdHMgKDMzKScsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICdCIGNlbGxzICgzNCknLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAnT3RoZXIgZW5kb3RoZWxpYWwgY2VsbHMgKDM1KScsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICdPdGhlciAoMzYpJywKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgJ090aGVyICgzNyknLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAnQiBjZWxscyAoMzgpJywKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgJ0VQQ0FNKyBjZWxscyBhbmQgY2hvbGFuZ2lvY3l0ZXMgKDM5KScpCiMgT3JkZXIgYnkgcm93IG5hbWUKYWxsQ2VsbHNNYXRyaXggPC0gYWxsQ2VsbHNNYXRyaXhbb3JkZXIocm93Lm5hbWVzKGFsbENlbGxzTWF0cml4KSksXQojIFNhbml0eSBjaGVjawpoZWFkKGFsbENlbGxzTWF0cml4KQpzcGVjaWVzRGF0YSA8LSAiYWl6YXJhbmkiCmBgYAoKQ29ycmVsYXRpb24gb2Ygd29vZGNodWNrIGxpdmVyIHdpdGggaHVtYW4gbGl2ZXIgZGF0YXNldCBmcm9tIEFuZHJld3MgKmV0IGFsLiogKDIwMjIpCgpgYGB7cn0KIyBMb2FkIGRhdGFzZXQKaHVtYW5OdWMgPC0gcmVhZFJEUygifi9Ecm9wYm94L1pvZS9zY2ZfdmVyc2lvbi9hbmFseXNpcy9jb3JyZWxhdGlvblRlc3RzL3NpbmdsZV9udWNfMjBfaHVtYW5fbWFwLnJkcyIpCiMgSXNvbGF0ZSBzaW5nbGUtbnVjIG9ubHkKaHVtYW5OdWMgPC0gVXBkYXRlU2V1cmF0T2JqZWN0KGh1bWFuTnVjKQpodW1hbk51YyA8LSBzdWJzZXQoaHVtYW5OdWMsIHN1YnNldCA9IGFzc2F5X3R5cGUgPT0gInNpbmdsZV9udWMiKQojIFJ1biBTQ1RyYW5zZm9ybQpodW1hbk51YyA8LSBTQ1RyYW5zZm9ybShodW1hbk51YykKIyBDaGFuZ2UgaWRlbnRzCklkZW50cyhodW1hbk51YykgPC0gaHVtYW5OdWNAbWV0YS5kYXRhJHN1Yl9hbm5vdGF0aW9uCiMgR2V0IGNsdXN0ZXIgYXZlcmFnZXMKaHVtYW5OdWNDbHVzdEF2ZXJhZ2VzIDwtIEF2ZXJhZ2VFeHByZXNzaW9uKGh1bWFuTnVjLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgYXNzYXlzID0gIlNDVCIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICBzbG90ID0gInNjYWxlLmRhdGEiKQpodW1hbk51Y0NsdXN0QXZlcmFnZXMkU0NUIDwtIG5hLm9taXQodChzY2FsZSh0KGFzLm1hdHJpeChodW1hbk51Y0NsdXN0QXZlcmFnZXMkU0NUKSkpKSkKIyBHcmFiIGdlbmUgbmFtZXMKbnVjR2VuZXMgPC0gcm93Lm5hbWVzKGh1bWFuTnVjQ2x1c3RBdmVyYWdlcyRTQ1QpCiMgTm93IHR1cm4gaW50byBsYXJnZSBkYXRhZnJhbWUKYWxsQ2VsbHNNYXRyaXggPC0gYXMuZGF0YS5mcmFtZShodW1hbk51Y0NsdXN0QXZlcmFnZXMkU0NUKQojIE9yZGVyIGJ5IHJvdyBuYW1lCmFsbENlbGxzTWF0cml4IDwtIGFsbENlbGxzTWF0cml4W29yZGVyKHJvdy5uYW1lcyhhbGxDZWxsc01hdHJpeCkpLF0KIyBTYW5pdHkgY2hlY2sKaGVhZChhbGxDZWxsc01hdHJpeCkKc3BlY2llc0RhdGEgPC0gImFuZHJld3MiCmBgYAoKCkNvcnJlbGF0aW9uIG9mIHdvb2RjaHVjayBsaXZlciB3aXRoIHdvb2RjaHVjayBQQk1Dcy4gRm9yIHRoaXMgY29ycmVsYXRpb24sIHJlYWQgaW4gdGhlIHdvb2RjaHVjayBsaXZlciBkYXRhc2V0IGF0IHRoZSBiZWdpbm5pbmcgb2YgdGhpcyBzY3JpcHQgYW5kIHRoZW4gcmVhZCBpbiB0aGUgd29vZGNodWNrIFBCTUNzIGJlbG93CgpgYGB7cn0KZ3JvdXBzIDwtICJnZW5lcmFsX2NlbGxfbGFiZWxzIgojIFN0YXJ0IHdpdGggbGl2ZXIgYW5kIHJlYWQgaW4gd29vZGNodWNrIFBCTUNzIGFnYWluCmxvYWQoIn4vRHJvcGJveC9ab2Uvc2NmX3ZlcnNpb24vYW5hbHlzaXMvaGVhbHRoeV9zYy9zZXVyYXRfb2JqZWN0cy9ub19kcm9wbGV0UUMvaW50ZWdyYXRlZF9QQk1DX2NjYV9rYW5jaG9yNV9zY0NsdXN0Vml6LlJEYXRhIikKSWRlbnRzKHNjU2V1cmF0KSA8LSAiaW50ZWdyYXRlZF9zbm5fcmVzLjAuNiIKIyBGaW5kIGNsdXN0ZXIgYXZlcmFnZXMKcGJtY0NsdXN0ZXJBdmVyYWdlcyA8LSBBdmVyYWdlRXhwcmVzc2lvbihzY1NldXJhdCwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICBhc3NheXMgPSAiU0NUIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICBzbG90ID0gInNjYWxlLmRhdGEiLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIGdyb3VwLmJ5ID0gZ3JvdXBzKQpwYm1jQ2x1c3RlckF2ZXJhZ2VzIDwtIGFzLmRhdGEuZnJhbWUobmEub21pdCh0KHNjYWxlKHQoYXMubWF0cml4KHBibWNDbHVzdGVyQXZlcmFnZXMkU0NUKSkpKSkpCiMgT3JkZXIgYnkgcm93IG5hbWUKYWxsQ2VsbHNNYXRyaXggPC0gcGJtY0NsdXN0ZXJBdmVyYWdlc1tvcmRlcihyb3cubmFtZXMocGJtY0NsdXN0ZXJBdmVyYWdlcykpLF0Kc3BlY2llc0RhdGEgPC0gIlBCTUMiCmBgYAoKIyMjIE91dHB1dCBwbG90cyBmb3IgYWxsIGNvcnJlbGF0aW9ucwoKYGBge3J9CiMgTm93IGZpbmQgaW50ZXJzZWN0aW5nIGdlbmVzCm1hdGNoZXMgPC0gaW50ZXJzZWN0KHJvdy5uYW1lcyhhbGxDZWxsc01hdHJpeCksCiAgICAgICAgICAgICAgICAgICAgIHJvdy5uYW1lcyh3b29kY2h1Y2tDbHVzdGVyQXZlcmFnZXMpKQojIExvb2sgYXQgaG93IG1hbnkgZ2VuZXMgbWF0Y2hlZApsZW5ndGgobWF0Y2hlcykKIyBNYWtlIG5ldyBtYXRyaWNlcyB3aXRoIG9ubHkgbWF0Y2hpbmcgZ2VuZSBuYW1lcwp0b0NvciA8LSBhbGxDZWxsc01hdHJpeFttYXRjaGVzLF0Kd29vZGNodWNrQXZlcmFnZXNDb3IgPC0gd29vZGNodWNrQ2x1c3RlckF2ZXJhZ2VzW21hdGNoZXMsXQojIERvIFBlYXJzb24KcGVhclZhbCA8LSBjb3IodG9Db3IsIHdvb2RjaHVja0F2ZXJhZ2VzQ29yLCBtZXRob2QgPSAicGVhcnNvbiIpCmhlYXRtYXAocGVhclZhbCkKICAgICAgICAjbWFpbiA9IHBhc3RlKCJQZWFyc29uIGNvcnJlbGF0aW9uIG9mIiwgc3BlY2llc0RhdGEsICJ2cyIsIHdvb2RjaHVja0RhdGEpLAogICAgICAgICN4bGFiID0gd29vZGNodWNrRGF0YSwKICAgICAgICAjeWxhYiA9IHNwZWNpZXNEYXRhKQogICAgICAgICNtYXJnaW5zID0gYyg2LDExKSkKI1Jvd3YgPSBOQSwKI0NvbHYgPSBOQSkKcGRmKHBhc3RlKCIuL2ZpZ3VyZXMvIiwgdGlzc3VlLCAiLyIsIHRpc3N1ZSwgIl8iLCBzcGVjaWVzRGF0YSwgIl9QZWFyc29uQ29yLnBkZiIsIHNlcCA9ICIiKSwKICAgIGhlaWdodCA9IDE1LCB3aWR0aCA9IDEzKQpoZWF0bWFwKHBlYXJWYWwsIG1hcmdpbnMgPSBjKDEzLDEzKSkKZGV2Lm9mZigpCiMgRG8gU3BlYXJtYW4Kc3BlYXJWYWwgPC0gY29yKHRvQ29yLCB3b29kY2h1Y2tBdmVyYWdlc0NvciwgbWV0aG9kID0gInNwZWFybWFuIikKaGVhdG1hcChzcGVhclZhbCkKICAgICAgICAjbWFpbiA9IHBhc3RlKCJTcGVhcm1hbiBjb3JyZWxhdGlvbiBvZiIsIHNwZWNpZXNEYXRhLCAidnMiLCB3b29kY2h1Y2tEYXRhKSwKICAgICAgICAjeGxhYiA9IHdvb2RjaHVja0RhdGEsCiAgICAgICAgI3lsYWIgPSBzcGVjaWVzRGF0YSkKICAgICAgICAjbWFyZ2lucyA9IGMoNiwxMSkpCiNSb3d2ID0gTkEsCiNDb2x2ID0gTkEpCnBkZihwYXN0ZSgiLi9maWd1cmVzLyIsIHRpc3N1ZSwgIi8iLCB0aXNzdWUsICJfIiwgc3BlY2llc0RhdGEsICJfU3BlYXJtYW5Db3IucGRmIiwgc2VwID0gIiIpLAogICAgaGVpZ2h0ID0gMTUsIHdpZHRoID0gMTMpCmhlYXRtYXAoc3BlYXJWYWwsIG1hcmdpbnMgPSBjKDEzLDEzKSkKZGV2Lm9mZigpCmBgYAoKCgoK